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Abstract 

O i 

The review of recently developed by the authors new techniques for covariant calculation 
of matrix elements in QED, the so-called formalism of "Diagonal Spin Basis" (DSB), is 
presented. In DSB spin 4-vectors of in- and out- fermions are expressed just in terms of 
their 4- momenta. In this approach the little Lorentz group, common for the initial and final 
states, is realized. This brings the spin operators of in- and out- particles to coincidence, 
allowing to separate in a covariant way the interactions with and without change of the spin 
states of the particles involved in the reaction and to follow in details the whole dynamics of 
the spin interactions. In contrast to methods of CALCUL group and others, the developed 
approach is valid for both massive fermions and massless ones. It is not necessary to introduce 
auxiliary vectors in DSB. Just 4-momenta of particles participating in reactions are required 



in it to construct the mathematical apparatus for calculations of matrix elements. We apply 
this formalism to the following processes: 1) Moller and Bhabha bremsstrahlung (e ± e~ — ► 
e e~7) in the ultrarelativistic (massless) limit when initial particles and photon are helicity 
polarized; 2) Compton back-scattering of photons of intensive circularly polarized laser wave 
focused on a beam of longitudinally polarized ultrarelativistic electrons (e + wjq — ► e + 7); 
3) e + e~-pair production by a hard photon in simultaneous collision with several laser beam 
photons (7 + n7o — ► e + + e~); 4) Bethe-Heitler process in the case of a linearly polarized 
photon emission by an electron with account for proton recoil and form factors; 5) the 
reaction ep — > ep'y with proton polarizability being taken into account in the kinematics when 
proton bremsstrahlung dominates; 6) orthopositronium 3-photon annihilation (e + e~ — > 3j). 
The results obtained with the help of the developed DSB-formalism certify its efficiency for 
calculating of multiparticle processes when polarization is to be taken into account. 

Introduction 

By now, the physics of spin phenomena has become an essential component of the research program 
at many large accelerators of new generation [1-5]. This is a consequence of, first, the successful 
development of the polarization technique, in particular, methods of obtaining polarized beams 
and advances in the construction of polarized targets and polarimeters [3] . Second, the electroweak 
interactions play an important role in the energy range of the current accelerators. They violate 
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both P and C parity [6,7], and also combined CP invariance [8]. This violation manifests itself in 
polarization effects, which are often used as precision tests of the Standard Model with accuracy 
unattainable in other experiments [1,5]. Third, it is necessary to go beyond the Standard Model 
to seek new particles and new types of interactions, and here spin can play a very important role. 

The advances in accelerator and polarization techniques have revealed new possibilities in 
the study of polarized particle interaction processes. Therefore, it is becoming more and more 
important to calculate theoretically the probabilities for various elementary-particle interaction 
processes, taking into account the particle polarizations and internal structure, and also to develop 
new computational tools. When the standard approach [9-12] is used to calculate the probabilities 
for various processes (i.e., to calculate the squared moduli of matrix elements), the inclusion of 
the particle polarizations greatly complicates both the calculations themselves and the structure 
of the expressions obtained. Their covariance is often lost. 

A natural way to simplify the calculations for reactions involving polarized particles is to 
calculate not the squared moduli of matrix elements, but instead the matrix elements themselves. 
This can be done in several different ways [13]. One way is to use the explicit form of the 
fundamental matrices and state functions written in a particular basis of the representation space 
of the Lorentz group in which they are defined. A noncovariant approach of this type has already 
been used for spin-1/2 particles by Powell [14] in 1949. The general theoretical development of 
this method is due to Sokolov [15]. This method continues to be used successfully to this day 
[16,17], owing to the appearance of powerful computer programs for analytic calculations. 

However, the most widely used method for calculating the matrix elements of QED processes 
is the covariant method not involving the use of the explicit matrices and wave functions. It was 
proposed in 1961 independently by Bellomo [18] and by Bogush and Fedorov [19]. This approach 
is based on the method of projection operators in elementary-particle theory developed by Fedorov 
[20]. 

The Bellomo method uses a trick which amounts to multiplying the matrix elements M 31 = 
^zQ^i of the transition from the initial state (\l/i) to the final state (^3), where Q is the interaction 
operator, by the quantity ^/ 1 Z^/ 3 /^/iZ^ 3 , so that the amplitude M 31 can be reduced to calculation 
of a trace :[] 

M 31 = Tr(P 31 Q) , P 3 i = *i ^3 , (1) 

P 31 = nZr 3 /(| *iZ# 3 I e*) , ^iZ^ 3 =| ^Z^ 3 I e** . (2) 

Here t\ and r 3 are the projection matrixs-diadics of the initial and final states [20]: r» = \Pj \Pi, (i = 
1,3). The operator Z in (2) is arbitrary. In Ref. 18 it was chosen to be: Z = 1. In recent years, 
the greatest progress in the development of the Bellomo method (in the ultrarelativistic, massless 
case) has been made by the CALCUL group [21]. The achievements of this group are widely 
recognized and extensively used by scientists all over the world. The CALCUL method has been 
generalized to fermions with nonzero mass in Refs. 22 and 23, but this generalization requires the 
introduction of additional vectors unrelated to the kinematics of the process under study, making 
it inconvenient to use. 

In the method proposed in Ref. 19, the operator P31 = ^1 \I/ 3 is constructed on the basis of a 
complex vector parametrization of the Lorentz group [24-26] and the operators T 3 i of representa- 
tions of this group in the space of particle wave functions [27,28] , which play the role of operators 

lr rhe indices have been chosen in accordance with future application of the results to the reaction 1 + 2 — ► 3 + 4 



for transitions from the initial to the final state: ^3 = T^^i , ^ 3 = ^iT 31 x . Here the operator 
P 31 = ^ l ~qj 3 is written as [13,19]: 

P31 = *i ^3 = riT^ 1 = T 3l l r 3 • (3) 



This version was originally developed for longitudinally polarized Dirac particles [29]. It was 
developed further by Fedorov [30-32] and his students (see Ref. 13 and references therein). In 
principle, the method developed by Fedorov (Ref. 13, Sec. 36) allows analytic expressions to be 
obtained for the matrix elements of various QED processes for arbitrarily polarized Dirac particles, 
either massive or massless, which is the main, decisive advantage of this method over that of the 
CALCUL group. However, the striving for generality is not always consistent with efficiency of 
the approach. 

For a number of QED problems, the development of the approach of Refs. 13 and 19 for 
calculating matrix elements of multiparticle processes is largely the result of progress in developing 
covariant methods of describing the spin properties of two-particle systems based on the use of a 
vector parametrization of the Lorentz little groups [13,37]. 

At the present time, the helicity basis introduced by Jacob and Wick [33] is very popular 
in high-energy physics. This is a consequence of the simplicity of the physical interpretation 
of helicity (the spin projection on the direction of the particle momentum), the fact that the 
center of mass of the system is distinguished in the helicity basis, and also the fact that the 
helicity amplitudes admit a simple partial- wave analysis using the 5*0(3) group [33]. In addition, 
studying the helicities of moving particles is analogous to studying the spins of particles at rest [13]. 
However, there are several important factors which prevent helicity from playing the dominant role 
in describing the spin projection of particles. One is that the helicity is not a particle characteristic 
which is invariant under Lorentz transformation [9,13]. Nevertheless, in the literature one can find 
articles with titles like "A Covariant Method for Calculating Helicity Amplitudes" (Ref. 34). In 
interpreting the dynamics of the spin interaction, amplitudes with and without change of the 
particle helicity are often referred to as amplitudes with and without spin flip. However, since 
the particle momentum is changed by the interaction, it is clear that such a classification is very 
arbitrary. Both types of amplitude actually describe a process with a change in the particle spin 
state. 

Many of these difficulties can be avoided for a particular choice of spin basis of a reaction, 
namely, the diagonal spin basis (DSB), in which the spin 4- vectors Si and s 3 of particles with 
4-momenta p\ and ^3 (sipi = s 3 p 3 = 0, s\ = S3 = —1) belong to the hyperplane formed by the 
4- vectors p\ and p 3 (Refs. 35 and 36): 

(v 1 v 3 )v 1 -v 3 (viv 3 )v 3 -v 1 

si = , , s 3 = - , (4) 

yWs) 2 - 1 ^{viv 3 y - 1 

where V\ = p\jvfi\ and v 3 = p 3 /m 3 . The spin vectors (4) obviously do not change under 
transformations of the Lorentz little group common to particles with 4-momenta p\ and p 3 [37]: 
Lpi, P3 Pi = Pi,L PltP3 p 3 = p 3 . We note that it will be a one-parameter subgroup of the rotation 
group with axis whose direction is determined by the vector [13,37]: 

o = c (Pi/Pio ~ Pz/Pm) , (5) 

where c is an arbitrary real number. The direction of a (5) possesses the property that the 
projections of the spins of both particles on it will have definite values even when the particles 



have different masses. Therefore, the DSB naturally makes it possible to describe the spin states of 
systems of any two particles (including ones with different masses) by means of the spin projections 
on the single common direction given by the vector (5).Q 

The fundamental fact that the Lorentz little group common to particles with momenta p% and 
£>3 is realized in the DSB leads to a number of remarkable consequences [35-42]. First, in this basis 
particles with 4-momenta pi (before the interaction) and p^, (after the interaction) have the same 
spin operators [38-40], which allows the covariant separation of the interactions with and without 
change of the spin states of the particles involving in the reaction, making it possible to trace the 
dynamics of the spin interaction. 

Second, in the DSB (4) the mathematical structure of the amplitudes is maximally simplified, 
owing to the coincidence of the particle spin operators, the separation of Wigner rotations from the 
amplitudes [35,36], and the decrease in the number of various scalar products of 4- vectors which 
characterize the reaction. Third, in the DSB the spin states of massless particless [p\ = p\ = 0) 
coincide up to a sign with the helicity states [40-42]. 

Use of the DSB does not lead to loss of generality, because the transformation to an arbitrary 
spin basis is carried out by means of Wigner D functions [43] . In the new expressions for the am- 
plitudes, the original amplitudes give the best representation of the dynamics of spin phenomena, 
and the D functions are purely kinematical in nature. 

Therefore, the DSB reveals new possibilities for developing methods to calculate matrix ele- 
ments and increasing the efficiency of such methods when the Bogush-Fedorov approach is used 
[13,19]. 

The first calculation of matrix elements in the DSB was performed in Refs. 35 and 36, using 
the spinor formalism. The amplitudes were calculated for the complete set of Dirac matrices 
Ti (i = 1,2 .. . 16) in which the arbitrary operator Q entering into (1) is expanded. We also note 
that the methods proposed in Refs. 18 and 19 are quite closely related, as was first shown in [44], 
where the various methods of calculating matrix elements were classified. 

Notation and abbreviations 

x = (x a ) is a three-dimensional vector, and x a (o = 1, 2, 3) are its components. 

P = (p k ) = {PoiP) is a four- dimensional vector in Minkowski space. 

xy = xiyi + X22/2 + £32/3 is the scalar product of the vectors x and y. 

P1P2 = Pifj,P2fj, = P10P20 — P~iP~2 is the scalar product of the 4- vectors p\ and p2- 

[xy\ is the vector product of the three-dimensional vectors x and y. 

e a bc is the three-dimensional Levi-Civita symbol. 

[xy\ a = e abc x b y c . 

(c)* b = e adb c d , c x x= [ex] , xc x = [x c] . 

x ■ y = (x a yb) is the dyadic formed from the vectors x and y . 

x ■ y = (x^y v ) is the dyadic formed from the 4- vectors x = (x^) and y = (y„). 

(x ■ y) = x ■ y + y ■ x , (x • y) = x ■ y + y ■ x are symmetrized dyadics. 

[x ■ y\ = x ■ y — y ■ x , [x ■ y}= x ■ y — y ■ x are alternating dyadics. 

(a x )^ = 1/2 e^a?" , aT = -a v » , (a = -a). 

2 The geometrical image of the difference of two vectors is the diagonal of a parallelogram, hence the name 
"diagonal spin basis" given by Fedorov. 



([a-b] x ) pu = e iiupa a p b u . 

[a, b, c] M = (\a-b] x c) p = e^ pa a u b p c a . 

g^v is the metric tensor in Minkowski space with signature (+, — , — , — ) 

^iivpa is the four-dimensional Levi-Civita symbol, e i23 = — 1- 

7 M are the Dirac matrices , a = a^ = (7a), 

,y5 __ — j _,0_,l_,2-,3 ~,5+ __ _,5^ 

The algebra of the Dirac matrices is: ab + ba = 2ab. 
abc = d — i~f 5 f, d = (ab + [a • b])c, f = ([a ■ b]) x c. 
For algebraic operations we use the notations: 
* for complex conjugation. 
+ for Hermitian conjugation. 
~ for the transpose, 
x for the dual 

for the dyadic product. 

QED stands for quantum electrodynamics. 
DSB stands for diagonal spin basis. 
OVB stands for orthonormal vector basis. 
RCS stands for real Compton scattering. 
VCS stands for virtual Compton scattering. 
CBS stands for Compton back-scattering. 

Everywhere we use the system of units in which the speed of light c and Planck's constant h are 
equal to unity: c — % — 1 . 

1 Spin operators in the DSB 

We shall use the following approaches in describing the spin properties of particles: (a) the ap- 
proach proposed by Bargmann and Wigner, in which the spin-projection operators are determined 
by using the generators of the Lorentz little groups. These are known in the literature as the 
Pauli-Bargmann-Lyubanskii operators [43,45]. (b) The covariant spin theory developed by Fe- 
dorov on the basis of vector parametrization of the Lorentz little groups and their representations 
[13]. These approaches are essentially equivalent. However, vector parametrization of the Lorentz 
group not only allows simplification of the theory of the spin properties of elementary particles, 
but also disposes of (see Ref. 13) commonly encountered, incorrect statements about some ap- 
proaches [9,10], such as "for a given momentum the spin projection on an arbitrary axis cannot 
have a definite value". 

We shall start from the fact that in momentum space the free state of a particle with 4- 
momentum p and spin projection 5 on the c axis is described by the state vector \p, 5 > (we drop 
the indices denoting the spin j, the mass m, and other particle characteristics). The particle spin 
j is defined as the angular momentum in the rest frame, where the orbital angular momentum is 
zero. It is therefore convenient to define the state vector \p, 5 > in terms of the state vector in the 
rest frame \p°, S >, where p° = (m, 0). Here we shall assume that the vector on which the spin is 
projected (i.e., the axis of spin projections c) in the particle rest frame is the spatial part of the 



spin 4-vector s° = (0, c), satisfying the conditions: s°p° = 0, s 02 = — c 2 = — 1. Let A p be a boost, 
i.e., a Lorentz transformation such that p = A p p°, s = A p s°, sp = 0, s 2 = — 1, s = (so, s), where 

c, s = — . (1.1; 



y m(po + m) J m 

Then 

\p,5>=T p \p°,5> , (1.2) 

where T p = T(A P ) is the operator for this transformation acting in the space of state vectors. The 
state vector \p, 5 > satisfies the equations: 

P»\p,6>=p»\p,5>, (1.3) 

a\p,5>=5\p,5> , (1.4) 

w 2 \p,5>=-j(j + l) \p,S> . (1.5) 

Here P M and o are the energy-momentum and spin-projection operators: 

a = s» uv > ( L6 ) 

where iu^ is the Pauli-Lyubanskii 4-vector [43]: 

«;„ = - ^ e^ v M"V , (1-7) 

and M^ A are the angular-momentum operators. Using (4), (1.6), and (1.7), we find that the 
spin-projection operators for the initial and final particles <J\ = ws\ and a 3 = WS3 in the DSB (4) 
coincide and have the form [38,39]: 

<n = a 3 = s^M^vlvl . (1.8) 

2^J(viv 3 ) 2 - 1 

It should be noted that in any other basis different from the diagonal one, the operators C\ and 
(T3 do not coincide and, therefore, do not commute with each other. 

The requirement that the Lorentz little groups coincide for particles with momenta p\ and p 3 
imposes rigorous constraints not only on the choice of particle spin vectors S\ and S3, but also on 
the spin-projection axes c.\ and c 3 (see (1.1)). As was shown in Refs. 35 and 36, c± and c 3 have 
the form: 

#31 _ #13 ( , n , 

I «31 I I «13 I 

where #13 (#31) is the spatial part of the relativistic difference of the 4-velocities of the first and 
third (third and first) particles (% = v i u,- = (%o,%)), defined as the velocity of the i-th 
particle in the rest frame of the j-th particle [36]: 

Vij = ViG Vj = A" 1 Vi . (1.10) 



Here A p x is the boost, A 1 i;j = v ° = (1,0), and i;^- = t>|j = 1, fjjo = ^jjo = ^j, and | iTy 



Vji l = y( v i v j) 2 ~ 1- The vectors #13 and -£7 3 i have the form [36]: 

«13 = «1 - U3 UlO - ^ ' U 31 = ^3 - V 1 \ V 30 - — . (1.11) 

To illustrate the properties of the DSB, let us consider an interaction process in the rest frames 
of the initial and final particles. In the first case [pi = (mi, 0)] the spin-projection axes C\ and c 3 
are parallel to the momentum of the final particle [this follows from (1.9) and (1.11)]: 

c 1 = c 3 = v 3 / \v 3 \ . (1.12) 

In the rest frame of the final particle [p 3 = (m 3 , 0)] the spin-projection axes are antiparallel to the 
momentum of the initial particle: 

ci = c 3 = -vxj I vi I ■ (1.13) 

Obviously, in these cases the Lorentz little group L PlP3 is a subgroup of the group of rotations 
about the direction of the momentum of the moving particle, which is the spin-projection axis for 
both particles. This is a special case of Eq. (5). 

Let us give another, equivalent representation for the spin-projection operator (1.6), expressed 
in terms of the antisymmetric matrix a(p) = [v ■ s] x , a(p) p — 0, and M^ u : 

a= l -{[vsY), u M^ . (1.14) 

In the DSB the alternating dyadics [v\ ■ si] and [v 3 ■ s 3 ] coincide: 

[vi ■ v 3 ] 



[vi ■ Si] = [v 3 ■ S3] 



y/(viv*) 2 ~ 1 



which ensures that the spin operators C\ and 0-3 coincide. We write the matrix a(p) 
expanded form: 

a(p) = - ( ? [&Pl ) (1 15) 

m \ -[cp\ Po(l - P ■ p/ ((po + m)Po)c) x J 

It is easily verifed that it has the same form in the rest frame (p = 0) and for c\\ p: 

«(P) = ( c- ) . ( L16 ) 

where c is an arbitrary unit vector in the first case and c = p/\p\ in the second. Therefore, study 
of the helicity states of moving particles is analogous to study of the spins of particles at rest, 
which is one reason for the popularity of the helicity basis. 

Let us now turn to spin-1/2 particles, the states of which are described by bispinors u 6 {p, s) 
satisfying the Dirac equation: 

(p — m)u (p, s) — 0, u (p, s)(p — m) — , (1-17) 



where u = u + j° with u s (p, s)u s (p, s) = m. The Dirac matrices satisfy commutation and recursion 
relations: 

7^ + 7^ = 2^, (1.18) 

YYY = g^Y + g up Y - g w i v + ii 5 e^ pa -? a , (1.19) 

yyy = g pv l 5 - i/2 e^ pCT 7 P 7 CT . (1.20) 
We write out these relations in the form without indices [13]: 

ah + ba = 2ab, abc — d — i^f > (1-21) 



1 5 ab = atrf 5 -i [a ■ b} x . (1.22) 

d = (ab + [a ■ b])c ,f = [a-b] x c=[a,b,c],f^ = e^a'tfc", (1.23) 



[a-6]x = l/2([a-6] x )^7^, (1-24) 

Let us also give some expressions which will be useful later on [13]: 



era — da = 2 cm = 2 (70:0), a x d — a a x = 2 a x a = 2 (^ya x a) , (1-25) 

where a — 1/2 ot llu Yl 1 ' -, and a is an arbitrary antisymmetric matrix. The first of these expressions 
can be obtained by multiplying (1.19) by l/2a fMU a p , then by l/2a fl a up , and subtracting the results. 
(The second one is found similarly.) 

In bispinor space the generators of the Loretz group M^ u have the form [43] : 

M^ = i/4 ( 7 V - 7'V)- (1.26) 

Then the spin-projection operator (1.14) for a spin-1/2 particle can be written as follows [32], 

using (1.26) and (1.22): 

i x 

a — — [v ■ s] x = - 7 5 s-0, aft = pa. (1-27) 

Therefore, the covariant electron spin-projection operator (1.14), which is directly related to the 
Loretz little group [13], differs by only the factor v from the widely used operator a (Refs. 9-12): 

a = l -Ys , [a'p]_ = , s 2 = -1 , sp = . (1.28) 

Here the commutation condition for the operator a and ft is satisfied automatically, as is easily 
verifed by using (1.25) and the equation [v ■ s] x p = 0. Therefore, for both moving particles and 
particles at rest the spin projection on an arbitrary axis can have a definite value [13]. The 
actions of the operators a (1.27) and a (1.28) on the particle state vector coincide because the 
Dirac equation is valid. 

Let us consider a binary reaction p\ + p 2 — > p% + J>4 in which particles 1 and 3 are of the same 
type, as are particles 2 and 4 (for example, electron-nucleon scattering ep — > ep, and so on). Given 
the spin structure of the matrix elements of this process, it is most convenient to use the DSB in 
which particles 1 and 3 and particles 2 and 4 have the same spin-projection operators. In order to 
construct the raising and lowering spin operators of the particles, we introduce the orthonormal 
vector basis (OVB) ua, with uaUb = gAB (A B — 0, 1, 2, 3) (Ref. 36): 

m = [n • n 3 ] x n 2 , n 2 = [p± ■ P3] x r/p , 



n 3 = ; = , n = , , (1.29) 

-(Ps-Pi) 2 J(Ps+Pi) 2 



where r is the 4- momentum of a particle participating in the reaction different from pi and p 3 , 
and p is determined from the normalization conditions n\ = n\ = n\ — —n\ — — 1. Therefore, the 
axes n and n 3 belong to the hyperplane formed by the 4- momenta p\ and p 3 , and n\ and n2 are 
orthogonal to them. The four vectors ha satisfy the relations [36]: 

[n A -n B ] x = 1/2 sab CD [n c • n D ] , [n A ,n B ,n c ] = -e A bc D n D . (1.30) 

They also satisfy the completeness relation: 

n ■ n - n x ■ n x - n 2 ■ n 2 - n 3 ■ n 3 = g , (1.31) 

by means of which an arbitrary 4-vector p can be written as: 

p = pn ■ no — pn\ ■ n\ — pn 2 ■ n 2 — pn 3 ■ n 3 . 

In the DSB, not only the spin-projection operators C\ and a 3 (1.27), but also the raising and 
lowering operators af and a 3 for particles 1 and 3 coincide. In the OVB (1.29) they have the 
form [39,40]: 

a = en = cr 3 = l/2 7 5 s 1 t) 1 = l/2 7 5 s 3 {) 3 = l/2 7 5 n n 3 = i/2fnh 2 , (1.32) 

a ±s = <jf s = af s = -l/2 7 5 n ±(5 , n ±5 = n 1 ± i5n 2 , 5 = ±1 , (1.33) 

au s {p t ) = 5/2u s ( Pl ) , a ±s u* s ( Pi ) = u ±s (p t ) , a ±s u ±s ( Pi ) = , (1.34) 

[aa ±s ]_ = ±5a ±s , [p l( r ±5 ]_ = [p 3 o ±& ]_ = , (1.35) 

where u 5 (pi) = u 5 (pi, Sj) are the bispinors of the first and third particles. 

Let us consider the projection operators for spin- 1/2 particles, r 6 = u s (pi) u s (pi) [13,20]: 

T s = l/4(p + m)(l-5^s) . (1.36) 

In the DSB the operators r/ (1.36) have the form [39,40]: 

4 = 1/4 (m + (£+n - £_n 3 ) + 5 7 5 (£_n - £+h 3 - m n 3 h )) , (1.37) 

r 3 5 = 1/4 (m + (£+n + £_n 3 ) - 5 7 5 (£_n + Z+h 3 + m n 3 n )) , (1.38) 



where £± = J(pip 3 ±m 2 )/2. Owing to (1.32), the spin parts of the projection operators for 
particles 1 and 3 can be made identical in the DSB, and so we have [46,47]: 

r/ = -1/8 fa + m) n s n* s , (1.39) 

where n* 6 = n\ — i5n 2 = n_$ and n$n} = —2. 

The bispinors of the initial and final states of the particles, u 5 (pi) and u s (p 3 ), can be related 
to each other by using the transition operators T 31 and T 13 = T 3 ^ [13,19]: 

u 5 (p 3 ) = T 31 u\ Pl ) ,u 5 (p 3 ) = u 5 ( Pl ) T 13 , (1.40) 



which in the DSB have the form [39,40]: 

1 + V 3 V 1 1 + ViV 3 



■ 31 



Ti 3 = , • (1.41) 



^/2( Vl v 3 + 1) ' '" ^/2( Vl v 3 + 1) 



Rewriting (1.41) in the OVB (1.29) and isolating the spin-projection operator a (1.32), we obtain 

[40]: 

T 31 = C~ CtV , T 13 = C + 2Cl 5 v , (1-42) 

from which we find the relation between the bispinors u s (p 3 ) and u s (pi) (Ref. 41): 

u S (p 3 ) = (C - 6^C) u\ Pl ) , u s ( Pl ) = (C + 5^C) u 6 (p 3 ) , (1.43) 

where £ ± = C,±/m. We also note that the Dirac equation can be used to reduce the transition 
operators T 31 and T 13 (1.41) to the same form [40]: 

T 31 = T 13 = h . (1.44) 

In the massless case the projection operators rf and t 3 (1.38) and (1.39) take the form [40-42]: 

t( =p 1 (l- 5 7 5 )/4 , t* = p 3 (1 + 5 7 5 )/4 . (1.45) 

It is easy to show that the operators rf and t 3 (1-45) satisfy the relations: 

^t* = 5 r* n 5 r! = -5 r* , (1.46) 

r( 7 5 = -5 4 , r 3 5 7 5 = 6 r 3 5 , (1.47) 

which imply that in the massless case the initial state is a helicity state, and the final state has 
negative helicity. 

Therefore, the DSB possesses a number of remarkable features which allow great simplification 
of the covariant calculation of the matrix elements for QED processes, to which we now turn. 

2 Calculation of matrix elements using the DSB 

The study of multiparticle reactions and the polarization phenomena arising in them requires 
effective computational tools. One is based on the use of the DSB (4). In the DSB the particle 
spin operators coincide. This allows the covariant separation of interactions with and without 
change of the spin states of the particles involved in the reaction. In the DSB, Wigner rotations 
[36,38], which are purely kinematical in nature, are separated from the amplitudes. This leads to 
maximal simplification of the mathematical structure of the diagonal amplitudes, and the resulting 
expressions give the truest reflection of the physical essential of spin phenomena. 

Let us turn to the calculation of the matrix elements of QED processes. They have the form 

M ±5 > s = u ±5 (p 3 )Qu 5 ( Pl ) , (2.1) 

where Q is the interaction operator, and u 6 {p\) and u ±5 (p 3 ) are the bispinors of the initial and 
final states, with u 5 (pi) u s (pi) = m , p\ = m 2 , (i — 1,3). 
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In the Bogush-Fedorov covariant approach [13,19] the calculation of matrix elements of the 
form (2.1) reduces to finding the trace: 

M ±s ' s = Tr(P± s > s Q) , p£ s < 5 = u s ( Pl ) u ±s (p 3 ) , (2.2) 

Psf = u 5 (pi) v*(p 3 ) = u 5 ( Pl ) u 5 ( Pl )T 13 = t*T 13 , (2.3) 

Psi' 5 = u 5 { Pl ) u ~ s (p 3 ) = a+ s u-\ Pl ) u ~ 5 (p 3 ) = a +& P^ & . (2.4) 

The operators P 31 ' determine the structure of the spin dependence of the matrix elements (2.1) 
in the case of transitions without spin flip M s,s and with spin flip M~ 6 ' 6 . Theyr explicit form in 
the DSB can easily be obtained by using (1.33), (1.37)-(1.41), and (1.44) (Refs. 39 and 40): 

^Pzi ={i+ + mh Q -^ n 3 n + tfy 5 ( £_ - m h 3 - £ + n 3 h )) , (2.5) 

AP 3 - L 5 ' S = -5 ( f- + m n 3 + £+ S 7 5 ) n s ■ (2.6) 

Equations (2.5) and (2.6) can be used to calculate the matrix elements, both with and without 
spin flip, for arbitrary Q. In particular, if the interaction operator reduces to the form 

Q = a + 7 5 6, (2.7) 

where a and b are 4- vectors, then for the matrix elements (2.1) we will have [39,40]: 

M 5 ' 5 = m (an + 5bn 3 ) , (2.8) 

M- 6 < 6 = -8£- (arts) + £+ (bn s ) . (2.9) 

Equations (2.5) and (2.6) can be written more compactly by using the operators (1.39) and (1.44), 
and also the expressions [39,40]: 

h 3 h h s = -5^ 5 ns , 7 5 n 5 n = 8h 3 n$ ,'y 5 n s n 3 = 5n ns . ( 2 -10) 

As a result, for the operators P 31 ' we will have [46,47]: 

4P 3 Y = (Pi + m) n s n h* 5 /2 , (2.11) 

AP 31 5 ' S = 5 (p! + m) h s n 3 . (2.12) 

Let us give yet another representation for the operators P 31 ' in (2.3), (2.4) in the DSB [42]: 

4Ps¥ = (Pi + m) ( 1 - ^ | (ft + m) , (2.5a) 

\V 2 (PiP3 + "i) \/2(p 1 p3-m 2 )y 

hd-W 5(pi + m) / 1 /, (Pi+Ps)^\ 

4^1 ' = < , r — m 



rn i I y2(p!p 3 - m 2 ) V PiPs + m 2 / 

<?7 5 ^ , (l>3-Pi)r \\,~ . x , OR x 

H , = [r + m } (p 3 + m) , (2.6a) 

^2( Pl p 3 + m 2 ) V PiPs - m 2 ; J 
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where rn\ is calculated by using the completeness relation (1.31): 

(rni) 2 = (rn ) 2 — (rn 3 ) 2 — r 2 , rn 2 = . 

Thus, the representation (2.5a) and (2.6a) is attractive in that it contains the Dirac operators 
only on terms of the particle 4-momenta pi,p 3 , and r, in contrast to (2.5), (2.6), (2.11) and (2.12), 
which involve h$ and fig. Moreover, the structure of the operators P 31 ' in (2.5a), (2.6a) is such 
that they automatically satisfy the Dirac equations: (p\ — m)P 31 ' = P 31 ' (p 3 — m) = 0. This was 
used to derive (2.5a) and (2.6a) from (2.5) and (2.6). 

Let us explain the choice of 4- vector r in terms of which the axes n x and n 2 (1-29) appearing in 
(2.5) and (2.6) are defined. First, it is chosen from the 4-momenta of the particles in the reaction 
under study, in contrast to the CALCUL approach, in which it is defined from considerations of 
convenience. Let us illustrate this for the example of the reaction e~(p\) + fi~(p 2 ) — > e~(p 3 ) + 
H~(p4) + 7(fc), which corresponds to Feynman graphs containing two fermion lines. For each of 
these lines it is necessary to construct the corresponding operators P 31 ' and P 42 ' in (2.5a), and 
(2.6a), expressed in terms of pi,p 3 ,ri and p 2 ,Pi,r 2 . For this process it is very convenient to make 
the choice: r\ = r 2 = k, so that kn 2 = kn 2 = and (kni) 2 = (ktio) 2 — (kn 3 ) 2 . The vectors r\ and 
r 2 can also be chosen to be the 4-momenta belonging to different fermion lines: r\ = p 2 , r 2 = p\ 
(here we are considering transitions p\ — > p 3 and p 2 — > P4). We note that the arbitrariness in the 
choice of the 4-vector r in (1.29), i.e., replacement of r by r', leads to the expression [42]: 

/ -t i iUi -e \ is<h r(ni ■ ni + iS[n ■ ns] x )r' 

\rri\)\r'u x ) 

and affects only the phase factor of the matrix elements M ±S ' S . 

Processes involving identical particles (for example, ee — > ee, ee — > ee-f, and so on) correspond 
to direct and exchange graphs [10]. They are associated with matrix elements M\ and M 2 of the 
form: 

Mi =u ± \p 3 )Q 1 u\ Vl ) -u ±s \p 4 )Q 2 u s \p 2 ) , 

M 2 = u ±5 ( P3 )Q3U 5 '(p 2 ) ■ u ±5 '(p,)Q,u 5 ( Pl ) , 

which are calculated as: 

Mi = TriP^Q,) Tr(P± 5 '> s 'Q 2 ) , M 2 = Tr(P^ 5 Q 3 P 4 ± /' 5 'g 4 ) ■ 

Therefore, the calculation of the direct graphs reduces to a product of traces, while that of the 
exchange graphs reduces to a trace extended by the product of the corresponding operators [13]. 
Let us give some useful expressions which are valid in the DSB [40]: 

a u 5 (pi) = (an + an 3 5^ 5 ) u s (p 3 ) + an s / y 5 W s (p 1 ) , 

a u s (p 3 ) = (an + an 3 5^ 5 ) u 5 (pi) + an s / y 5 u~ s (p 3 ) , 

u*(pi) a = u s (p 3 )(an - an 3 5rf) - an* s u~ 5 {pi) 7 5 , (2.13) 

u 5 {p3) a = u 5 (pi)(an - an 3 5rf) - an* s u~ 5 (p3) 7 5 , 
where a is an arbitrary 4-vector {n* s = n_$). 
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In the DSB (4) the particle spin vectors are expressed in terms of the 4-momenta, so that the 
number of independent scalar products entering into the final expressions for the matrix elements 
after calculation of the traces (2.2) is decreased. For the same reason the circular-polarization 
vector e\ of a photon with 4-momentum k emitted by a particle in the transition p\ — > p 3 is 
conveniently defined by using the 4- vectors Pi,p 3 and k (Refs. 21 and 40): 

[n ■ n 3 }k + i\ [n ■ n 3 ] x k [Pi ■ Ps] , n , A , 

«a = Wp , [«o ■ »,] = -^ , (2.14) 



where p — J— ([pi • p 3 ]k) 2 /2t; + £-. Then for the dyadic e\ ■ e* x we easily find: 

1 ( k-nx+ni-k k-k [£:-ni] x \ 

eA -^ = 2 \- g+ — ^ — + yq +tX ^T) • (2 - 15) 

Using (1.25) and (1.22) the operators e±\ (e* x = e^\) can be written as follows [40]: 

e±A = N 13 (kpsPi(l T A 7 5 ) - PaPiMl ± A 7 5 ) =F 2p i p 3 X-y 5 k) , (2.16) 

N^ 1 = 2 1 / 2 (8 PlP3 ■ Pl k ■ p 3 k - m 2 {{2p l kf + {2p 3 k) 2 )f' 2 . 
In the massless case (p 2 — p\ — 0) the operators P 31 ' in (2.5) and (2.6) take the form [40]: 

AP 5 3 f = £(1 + o 7 5 )(l + n h 3 ) , AP^ 5 = -^(1 + 5 7 5 )n 5 , (2.17) 



where £ = £ + = £_ = JpiPz/2. Similarly, from (2.5a) and (2.6a) we have: 

API? = i 1 + */)*& , AP 31 S > S = -^ (1 ~ Wl ^ 3 . (2.17a) 

Using (1.25), it is easy to show that the representations (2.17) and (2.17a) are equivalent. As noted 
above, in calculating processes with the emissions or absorption of a real photon with 4-momentum 
k, it is convenient to make the choice r = k for the 4- vector r entering into (2.17a). Then the 
denominator of the operator P 31 ' in (2.17a) takes the form m\^/2p\p 3 = \flp\k ■ 2p 3 k, and we 
obtain a result similar to that of Ref. 23, except that our expressions involve the 4-momentum of 
a real photon, and not an auxiliary lightlike 4-vector k. These points are very important for our 
approach, in which we use only the 4-momenta of the particles participating in the reaction. 
Using (1.46) and (1.47), Eq.(2.13) can be written as [40]: 

a u (pi) = {ario — an 3 ) u (p 3 ) — 5arisu~ (pi) , 

a u6 (Ps) = ( an o + an 3 ) u 5 (pi) + 5arisu~ s (p 3 ) , 

vf(pi) a = (an — an 3 ) m 5 ^) — 5an* s u~ s (pi) , (2-18) 

u?{P3) a = ( a ^o + a n 3 ) w 5 (pi) + San* s u~ 5 (p 3 ) . 

In the massless case the relation between the bispinors of the initial and final states takes a 
particularly simple form [see (1.44)]: 

Ap3) = Tib " S(Pl) ' Api) = Jk Ap3) ■ (2 ' 19) 
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In this massless limit, the terms containing 7 5 /c in (2.16) can be dropped, owing to gauge invariance. 
As a result, for the operators e±\ we obtain the expressions used by the CALCUL group [21,48]: 

e ±A = N 13 (kp 3 Pi(l T A 7 5 ) - PsPikO- ± A 7 5 )) , (2.20) 

N^ 1 = 4 (pip 3 -pik-p 3 k) 1/2 . 
Using (2.18)-(2.20), we can easily verify the correctness of the expressions [21,40,48]: 

ex u\ Pl ) = -(1 + SX) 2 Pl k N 13 p 3 u\ Pl ) , 

rf(p 3 ) e* x = (1 - SX) 2p 3 k N 13 tf(p 3 ) px . (2.21) 

If photon emission occur in transition pa — >■ p B , then, making the replacement (pi, p 3 ) — ► (pa,Pb) 
in (2.20), we obtain the operators <3aab, whose action on bispinors is the same as that of <3a = e A13 
except for a phase [21,40]: 

<5ai3 = £\ab exp (i(pAB) , exp {i^ab) = iX2 ll2 e xl3 n 2 ( A B) , (2.22) 

where U2(ab) are unit vectors: 

n 2 (AB) = [PA-PB] X k/p (A B) , P(AB) = (2p A PB • PAk ■ p B k) 1/2 . 

Up to now our discussion has pertained to the case with only electrons in the initial and final 
states. If one state is an electron and the other a positron, the amplitude of the process will have 
the form [9]: 

M 31 *' s = { ^ { r P d Q ^ (Pl \ (2-23) 

3 I u ±s (p 3 )Qu 5 (-pi) , v ; 

where ■u <5 (— p\) and u ±5 {— p 3 ) are the positron bispinors in the final and initial states, with 
u 6 \—pi)u 6 (—pi) = —m (i = 1,3). The upper amplitude in (2.23) corresponds to pair annihi- 
lation, and the lower one to pair production. To construct the operators 

^f' 5 = u\ Pl ) u ±5 (-p 3 ) , P£ s ' 5 = u\- Pl ) u ±5 (p 3 ) , (2.24) 

used to reduce the determination of the matrix elements (2.23) to calculation of the traces M ±s,s = 
Tr(P 31 ' Q), we need to use the relation between the positron and electron in the DSB [36,38]: 

u\-p) = -5 7 5 u-\p) , u 6 (-p) = u-\p) 5 7 5 . (2.25) 

As a result, for the operators P 31 ' used to calculate the amplitudes for pair annihilation, we 
obtain: 

AP 3 [ = S (pi + m) n hs , AP 3l ' = - (pi + m) n s n 3 h* s /2 . (2.26) 

Similar expressions can be obtained for the operators P 31 ' in the case of pair production. 

We have used this formalism for calculating matrix elements in the DSB to obtain the cross 
sections for several real QED processes, to which we now turn. 
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3 The cross sections for the processes e ± e — » e ± e 7 
in the ultrarelativistic massless case 

Moller and Bhabha bremsstrahlung e ± e~ — » e ± e _ 7 are background processes in studying hadron 
states. Moreover, the study of these processes allows verification of QED in higher orders of per- 
turbation theory. The cross sections for these processes are quite awkward, even in the ultrarela- 
tivistic limit. Only relatively recently has it been possible to write them down in a compact form 
for unpolarized [49] and transversely polarized initial particles [50] . Using the methods described 
above [Eqs. (2.17)-(2.22)], the present authors have obtained [40] compact expressions, in the 
ultrarelativistic, massless limit, for the differential cross sections of the processes e ± e~ — > e ± e~ r y 
for the case where not only the initial particles but also the photon are helically polarized. As 
was shown in Ref. 40, the cross sections for these processes are written as the product of two 
factors, one universal and coinciding with that obtained earlier [49] for unpolarized particles. Let 
us consider Moller and Bhabha bremsstrahlung, 

e"(pi) + e ± (p 2 ) -> e-(p 3 ) + e ± (p 4 ) + l(k) , (3.1) 

assuming that the initial and final e ± particles are massless (p 2 = 0, % = 1,2,3,4). The details of 
the calculations of the matrix elements for (3.1), which correspond to eight Feynman diagrams 
[9,10], are given in Ref. 40, and so we shall not dwell on them here. We introduce the invariant 
variables [49,50]: 

s = (pi+ p 2 ) 2 ,t = (pi- p 3 ) 2 , u = (pi - p A f , 

s = (p 3 + Pa) 2 , t = (p2 - Pa) 2 , u = (p 2 - Pzf , (3.2) 

and also the notation 5, 5 and A for the helicities of the initial particles and the photon, respec- 
tively. Then the differential cross sections for the processes e~e ± — > e~e ± 7 in the case helically 
polarized initial leptons and photon have the form [40]: 

o; 3 a 3 

dcr M = —z- A M W M dT , da B = — - A B W B dT , (3.3) 

tt z s n s 

A M = A MB /t t uu , A B = A MB /t t s s , (3.4) 

A MB = 1/2 {ss'(s 2 + s' 2 ) + tt{f + t' 2 ) + uu{u 2 + u' 2 ) 

+ 5 5' (ss'(s 2 + s' 2 ) - tt'(t 2 + t 2 ) - uu'(u 2 + u 2 )) (3.5) 

+ 5 \ (ss'(s 2 - s' 2 ) - tt(t 2 - t 2 ) - uu'(u 2 - u' 2 )) 
+ S'X (ss'(s 2 - s' 2 ) + tt{t 2 - t 2 ) + uu'(u 2 - u' 2 ))} , 

^" = -(m + m"m~^J ■ (3 - 6) 

\P\k p 4 k p 3 k p 2 k) 

4 j,d 3 p 3 d 3 p A d 3 k 

dT = 5 (p 1 +p 2 -p 3 -p 4 - k)- — - — — , 

2p 3 o 2p 40 2u 
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where xa = PAk , (A = 1, 2, 3, 4), and a is the fine-structure constant. We note that Eq. (3.5) for 
Amb is invariant under crossing transformations: 

Vi ^ —Pi , x 2 <-> — %i , s ^ u , s <-> u ,5 ^ —5 . (3.8) 

The expressions for Amb, Wm, and W B can be written in a different form [40]: 

A MB = 1/2 {(1 + 55') ((1 + 5\)ss's' 2 + (1 - 5\)ss's 2 ) + 

+ (1 - 55') ((1 + 5'X) {tit 2 + uu'u 2 ) + (1 - 5' A) (tt't' 2 + uu'u' 2 ))} , (3.9) 

-W M = ^ + ^ + ^ + ^ + ^L_ + JL_^ (3.10) 

X1X2 X3X4 X1X3 X2X4 X1X4 ^2^3 

s s t t U U /„ ,,n 

W B = + + + . (3.11) 

X1X2 X3X4 X1X3 X2X4 X1X4 x 2 x 3 

When soft photons are emitted (s = s , t = t , u = u), A M and A B take the form [40]: 
(5 = ±1,5' = ±l) 

s 2 + u 2 s 2 + t 2 2s 2 ,,,{s 2 -u 2 s 2 -t 2 2s 2 \ .. 

^ = ^^ + ^^ + ^ + M (^^ + ^^ + ^J ' (3 ' 12) 

u 2 + s 2 u 2 + t 2 2u 2 ^,(u 2 -s 2 u 2 + t 2 2u 2 \ .„.,„, 

^ = ^^ + ^^ + lr -«(^^ + ^^ + lr ) • (3.13) 

They differ only by overall factors from the cross sections for elastic processes e ± e~ — > e ± e" when 
the initial particles are longitudinally polarized (see Ref. 10). 
For unpolarized photons, from (3.9) we have: 

A MB = (1 + 55') ss (s 2 + s' 2 ) + (1 - 55') (ti(t 2 + i 2 ) + uu'(u 2 + u' 2 )) . (3.14) 

Therefore, the ratio of the cross sections for particles with parallel and antiparallel spins have 
the same form [40] for the two reactions e ± e~ — > e ± e~7 (as in the case of the elastic processes 
e ± e~ — > e ± e~; see Ref. 10): 

da v _ tt'jt 2 + t' 2 ) + uu'ju 2 + u' 2 ) 
da-i[ ss' (s 2 + s' 2 ) 

4 Polarization phenomena in the three-photon annihila- 
tion of orthopositronium 

In recent years the three-photon annihilation of orthopositronium 3 S\ — ► 37 has attracted a great 
deal of attention, because experiments to measure the decay width of orthopositronium revealed 
a discrepancy with the theoretical predictions [51]. Several attempts have been made to resolve 
this contradiction. In Ref. 52, relativistic corrections were included in the cross section for the 
annihilation of a slow e + e~ pair into two or three photons, but this did not solve the problem. The 
contribution of the five-photon decay mode of orthopositronium, calculated in Ref. 53, indicates 
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that this mechanism also cannot eliminate the discrepancy in the width. All these problems, 
including the results of Refs. 51 and 53, require further analysis and confirmation. The work of 
Ref. 54 does not represent an attempt to resolve the orthopositronium problem. Almost all the 
known results pertaining to polarization phenomena in the reaction e + e~ — »■ 37 were obtained 
there, but by calculating the matrix elements in the DSB. The purpose of that study was to 
demonstrate the effectiveness of that method for a process to which the CALCUL method is 
inapplicable. A key feature of the technique is the very specific choice of the photon polarization 
vectors (2.20), which is valid only for the massless case. 

The main process determining the positronium lifetime is three-photon annihilation. Here the 
decay probability can be related to the cross section for annihilation of a free pair [10]: 

e"(pi) + e + (p 3 ) -► 7 (fci) + 7 (fc 2 ) + l(fa) • (4.1) 

Since the momenta of the electron and positron in positronium are small [10] (\ Pi |, | p 3 | ~ ma, 
where a is the fine-structure constant), in calculating the annihilation cross section they can be 
considered to be at rest at the origin [i.e., we assume that Pi = P3 = p = (m, 0, 0, 0)]. In this case 
the matrix element of the reaction (4.1) takes the form: 

M^ s = u ±5 (-p)Qu 5 (p), (4.2) 

where u 5 (±p) are the electron and positron bispinors, u & (±p)u s (±p) = ±m, and Q is the interac- 
tion operator, which corresponds to six Feynman diagrams [10]. Let us consider the kinematics 
of the process e + e~ — > 3 7 in the e + e~ cm. frame, in which the momenta p\ and p 3 have the form 
P\ = (po, 0, 0, —ma) and p 3 = (po, 0, 0, ma), po = my/1 + a 2 . We introduce the OVB a^- 

a = (1, 0, 0, 0) , 0l = (0, 1, 0, 0), a 2 = (0, 0, 1, 0), a 3 = (0, 0, 0, 1) , 

using which we find : 

pi = -£_a 3 + £+a , P3 = £-03 + £+a , (4.3) 

Si = £+a 3 - £-a , s 3 = £+a 3 + C^o , (4.4) 

where £ ± = £±/ra, £ + = yl + a 2 , and £_ = a, with sipi = s 3 p 3 = 0, and sf = s 2 = — 1. 
Therefore, in the limit a — > the electron and positron spin vectors S\ and s 3 (4.4) in the DSB 
(4) coincide: 

si = s 3 = a 3 , (4.5) 

i.e., the direction of the positron motion is singled out as the common axis of spin projection. The 
momentum conservation law 

ki + fa + h = (4.6) 

determines the annihilation plane in which the photon momenta lie. We shall also assume that the 
vectors ai and a 3 lie in this plane, while the vector a^ is normal to it, i.e., a^n^ = , Hi — ki/uii, 
and Hi 2 — 1, (i — 1, 2, 3). 
Let us construct the photon circular-polarization vectors e\i = (0,e\i): 



e\i = ([a 2 Hi] + iXi a 2 )/V2 , e Xi Hi = , e Ai e A * = 1 , (4.7) 

where Aj are the photon helicities, A« = ±1. 
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In the limiting case that we are considering, the operators (2.26) used to calculate the matrix 
elements (4.2) have the form: 

4P 3 Y = 5 (m + p) a s , 4P 3 ^' 5 = - (m + p) a 5 a 3 a* 5 /2 , (4.8) 

where a±s — Oi ± i5a 2 , 5 = ± 1 (a} = as). 

The explicit form of the matrix elements M ±S ' S for the process (4.1) in the case of circularly 
polarized photons was obtained in Ref. 54: 

3 

2 3/2 M 5,5 = 6J2<Xi(S\ i + c i ) (n jk - l)/ro , (4.9) 

i=l 



2 3/2 M -«5,«5 = Y t a i s i (n jk - l)/m , (4.10) 

where ctj (i=l,2,3) are the polarization factors: 

ai = (1 + A 2 A 3 )(1 - AiA 2 ), a 2 = (1 + AiA 3 )(l - A 2 A 3 ), a 3 = (1 + AiA 2 )(l - A1A3), (4.11) 

«i «2 = «l «3 = «2 «3 = j a^ = 4 q, , (4-12) 

and the quantities q, s«, n^ = ^fcj (i,j, k — 1, 2, 3 ) are: rijfc = WjWfc = CjCfc + SjS*;, Q = <J 3 nj, s« = 
aiflj, s 2 + c 2 = 1, with the indices i, j, and fc in (4.9) and (4.10) representing a cyclic permutation 
of the numbers 1,2,3. 

The matrix elements (4.9) and (4.10) determine the annihilation of a free e + e _ pair in the case 
of parallel (M s,s ) and antiparallel (M~ s,s ) spins of the electron and positron. They are real and 
vanish if all the photons have the same helicity, i.e., when Ai = A 2 = A3. 

The differential cross section for the process (4.1) is expressed in terms of the matrix elements 
M ±5 ' s in (4.9) and (4.10) as 

^ = §^(4 #)**<*- *.-*■-*»). («3) 

where v is the relative velocity of the e + and e~ in the cm. frame (v ~ a). We introduce the 
notation (T Xl ^ e l M = 1/2 ( (M 5 > 5 ) 2 + (M 5 >~ 5 ) 2 ). Then for <J Xl £l M we find [54]: 



a 



Ai,A2,A3 

i=l 



2 a, ( 1 + 5A, Ci ) (1 - n jk f/2m 2 , (4.14) 



which determines the annihilation cross section when all the particles except the electron are 
helically polarized. For the quantities 

a \iMM = i/4]T((M 5 ' 5 ) 2 + (M 5 -- 5 ) 2 ), a£ + = 1/2 J2 i(M 5 ' 6 ) 2 + (M* 5 '" 5 ) 2 ) , 

8 A2,A3 

the meaning of which is clear from the notation, we find 

<7A1.A2.As = ( ai (1 _ n23 )2 + « 2 (1 _ ni3 )2 + « 3 (1 _ n 12 ) 2 )/2m 2 , (4.15) 
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(r£ + = 2 ((1 + a lCl )(l - n 23 ) 2 + (1 - a lC2 )(l - n 13 ) 2 + (1 - a lC3 )(l - n l2 f)/m 2 . (4.16) 
In the case of unpolarized particles we obtain the well known result [10] 

w= £ ^1.^.^3 = 4 ( (1 - m 2 ) 2 + (1 - nis) 2 + (1 - nas) 2 )/m 2 . (4.17) 

Let us calculate the probability for the process (4.1) when one of the photons is linearly polar- 
ized in the annihilation plane (a x ) or perpendicular to it (a y ) (and the other two are unpolarized), 
and also the degree of linear polarization pi\ 

Pi = (cr y ~ cr x )/(a y + a x ) . (4.18) 

For this we go from the helicity states | +1 > and | — 1 > of a photon of momentum k\ 

| +1 >= (| x > +i | y >)/v / 2 , \ -1 >= (\ x > -i \ y >)/v / 2 , 
to states with linear polarization | x > and | y >: 

< x |= (< +1 | + < -1 \)/V2 , <y\=i(<+l\-<-l \)/V2 . 
Then for the amplitudes and probabilities we find [54] 

M 5 X = (M A2 5 ;/ 3 + M_\ 2 S X3 )/V2 , M-J = (M^J + MI* *« 
M S y = i (mJ;/ 3 - M_%%)/y/2 , M- 5 = i (M^ - MZ^ 



a x = 1/4 Y, (I Mi | 2 + | M~ 5 | 2 ), a y = 1/4 £ (| M s y | 2 + | M y ~ 5 | 2 ) . 
5A2A3 (5A2A3 

Computing cr x , a y , and pi, we obtain [54] 

a x = 2/m 2 (A-B), a y = 2/m 2 (A + B), Pl = B/A , (4.19) 

A = (1 - n 12 ) 2 + (1 - n 13 ) 2 + (1 - n 23 ) 2 , B = (1 - n 12 )(l - n 13 )(l - n 23 ) . 

Equation (4.19) coincides with the result of Refs. 55 and 56. 

Let us use (4.9) and (4.10) to construct the amplitudes for orthopositronium annihilation [10]: 
Xl,i = M++,X li0 = (M+~ + M~ + )/y/2, and X 1; _! = M~(X li±1 = X M ), corresponding to the 
projections of the total spin of the system on the direction of a 3 equal to +1, 0, and -1, and the same 
for parapositronium (with total spin and projection equal to zero): X 0i o = (M +_ — M _+ )/ ' \[2. 
We find [54]: 

X liS = M 5 ' 5 , X lfi = s/2 M 5 >- s , X ,o = . (4.20) 

Summation of X 2 S and X 2 over the photon polarizations A 2 and A 3 gives 



Y, X b=-!k ( ( ' + rt 'Vi) 2 d - "2,r + CI - 5AiC 2 ) 2 (] „ 1:1 ) 2 + (I ajrO-Cl „ 12 ) 2 ). 

A2,A3 



m? 



£ X 2 = 4 ((1 - n 12 ) 2 (l - c 2 ) + (1 - n 13 ) 2 (l - c 2 ) + (1 - n 23 ) 2 (l - c\))/m 2 . 

^2,A3 
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Averaging the squares, we again obtain the well known result [10] 

£ (X 2 1>+1 + X% + X$ t _ x ) =Aa, 

Ai,A2,A3 

where a is given by (4.17). 

It was shown in Ref. 55 that the amplitudes for the three-photon annihilation of orthopositro- 
nium Hfi (Xl j+ i, Xl i0 , Xl,-i) can be written as 

Hfi —tu , u = U1 + U2 + U3 , (4-21) 

where the vector u\ is a function of the photon polarization vectors: 

u x = ei (e 2 e 3 - e 2 'e 3 ') + e/ (e 2 e 3 ' + e 3 e 2 ') , (4.22) 

e, = e Xi , e/ = [e^] (« = 1, 2, 3) , 

and the vectors u 2 and 7/3 are obtained from u\ by cyclic permutation of the indices. The complex 
vector t characterizes the triplet state of orthopositronium. 

Let us construct the tensor $ = u ■ u * in terms of which the three-photon annihilation prob- 
ability is expressed. According to our calculations, the tensor $ can be written as three terms 
[54]: 

$ = (H&! + a 2 $ 2 + a 3 $3 , Tr($) = 01 + a 2 + a 3 = m 2 a XlM ^ 3 /2 , (4.23) 

01 = «i(l - n 23 ) 2 , a 2 = a 2 (l - n i3 ) 2 , a 3 = a 3 (l - n i2 ) 2 , (4.24) 

where each of the tensors $j (i=l,2,3) is just the "beam tensor" (the three-dimensionally covariant 
polarization density matrix) of the corresponding circularly polarized photon [57]: 

$i = e Xl ■ e A * = (1 - Hi ■ n t )/2 + i/2Xi n t x , $,n, = , Tr^) = 1 . (4.25) 

Since the tensor $ corresponds to the sum of three waves and its trace Tr(§) = a,\ + a 2 + a 3 
coincides up to an overall coefficient with the probability o" Al ' A2 ' A3 , each of the Oj (i = 1,2,3) in 
(4.24) determines the probability for the appearance of a single photon having polarization vector 
e\i and direction of motion Hi. 

5 The reaction ep — » epj and the proton polarizability 

There has recently been much interest in studying Compton scattering on nucleons at low and 
intermediate energies. This is because the fundamental structure constant of the nucleon-the elec- 
tric and magnetic polarizabilities- can be determined in this process. The nucleon polarizabilities 
contain important information about the nucleon structure at large and intermediate distances, 
in particular, the radius of the quark core, the meson cloud, and so on (see the detailed discussion 
of these questions in Refs. 58 and 59). Knowledge of the amplitudes for Compton scattering 
on nucleons is also required to interpret the data on photon scattering on nuclei. For example, 
studies of this type can answer the question of how greatly the electromagnetic properties of free 
and bound nucleons differ. 

All the experimental results on the proton polarizabilities have been obtained from data on 
elastic '-/p scattering below the pion photoproduction threshold. However, it has recently been 
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shown that measurement of the proton polarizabilities at the Novosibirsk storage ring with electron 
beam energy of 200 MeV using an internal jet target appears very promising. As proposed in [60], 
this can be done using the reaction 

e"(pi) +p + (qi) -> e~(j9 2 ) +P + (<?2) + l{k) (5.1) 

in the kinematics corresponding to electron scattering at small angles and photon scattering at 
rather large angles, which corresponds to small 4-momentum transfer from the initial electron to 
the final 7 and p. In the lowest order of perturbation theory, the process (5.1) is described by the 
three graphs shown in Fig. 1. 

k> n k 

P2 Pi < P2 
► 1 - 1 — ► 



q 2 q! — q 2 




Figure 1: Graphs corresponding to the reaction ep — > epy. 

The first two (a) and (b) correspond to electron bremsstrahlung (Bethe- Heftier graphs), and 
the third (c) corresponds to proton bremsstrahlung [graph with virtual Compton scattering (VCS) 
on a proton]. The kinematics described above was chosen for the following reasons. First, the 
subprocess of real Compton scattering (RCS) on the proton is realized in it, because at small 
electron scattering angles the virtual photon with 4-momentum r = p\ — P2 (see Fig. 1) becomes 



almost real. Here the quantity |r| = J—(pi — P2) 2 turns out to be small, \r\ ~ m, where m is 
the electron mass. Second, for electron scattering at small angles and photon scattering at fairly 
large angles, the contribution of the graph corresponding to proton bremsstrahlung dominates, 
i.e., it is several orders of magnitude larger than the contribution of the Bethe-Heitler graphs to 
the cross section for the process (5.1) (Ref. 61). This is the main requirement needed to isolate 
the subprocess of Compton scattering on a proton [60] from the reaction ep — > epj. 

The estimates made in Ref. 60 using the method of equivalent photons and the scalar model 
showed that the reaction (5.1) offers a good possibility of obtaining high-statistics data on the 
Compton scattering cross section and the proton polarizability. Measurement of the electric and 
magnetic polarizabilities of the proton (a p and /3 p ) with higher accuracy than in earlier studies is 
one of the most important problems to be solved by experiments in the near future [62,63]. 

However, to obtain high-statistics data on the cross section for 7p scattering and the proton 
polarizability it is essential to use a theoretical model more accurate than that in Ref. 60. It must 
include both the spin properties of the particles and the main structural parameters characterizing 
the electromagnetic structure of the hadron. The model can be based on the result of Ref. 64, 
where a general calculation of the reaction ep — > ep~y was performed. The cross section was 
expressed in terms of 12 form factors corresponding to the VCS subprocess on the proton (i.e., 
the contribution of the graph in Fig. lc) and two form factors corresponding to the Bethe-Heitler 
graphs. 
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The differential cross section for the reaction ep — > ejrj in the above kinematics was calculated 
in Ref. 65. It was expressed in terms of the six invariant amplitudes for RCS [58,66], and also the 
electric and magnetic form factors of the proton [10]. 

Let us consider the amplitudes corresponding to the graphs of Fig. 1. The sum of the two 
Bcthc-Heitler graphs (a) and (b) corresponds to the matrix element 

Mx = u( P2 )Q>(pi) ■ u(q2)T II (q 2 )u(q 1 ) 1 , (5.2) 

q 2 

n „ fl Pi-k + mp 2 + k + m 
—2pik 2p 2 k 

TM 2 )=hl, + -^h{qi,-l,q), (5.4) 

where u(pi) and u(qi) are the bispinors of electrons and protons with 4-momenta pi and qi, p 2 = 
m 2 , q 2 = M 2 , u(pi) u(pi) = 2m, u(qi) w(%) = 2M, (i = 1,2), /j, p , f ± , and f 2 are respectively the 
anomalous magnetic moment and the Dirac and Pauli form factors of the proton [10], q = q 2 — q\ is 
the momentum transfer, e is the polarization 4- vector of a photon with momentum k, ek = k 2 = 0, 
and M is the proton mass. 

In the limit of interest |r| ~ m, the matrix element corresponding to the graph of Fig. lc will 
be expressed in terms of the six invariant RCS amplitudes Tj (i = 1,2. . .6) obtained from the 
theory of dispersion relations and the data on 7r-meson photoproduction on nucleons [66] . It has 
the form [64] 

M 2 = u(P2)Yu( Pi ) ■ u(q 2 )M IIU e u u(q 1 ) — , (5.5) 



M^ = ^ (T\ + T 2 K) + ^~— (T 3 + T,K) + 



C,C U (T+Tk) , D^Dv 



The tensor M^ v is constructed using a set of four mutually orthogonal 4- vectors C, D, B, and K: 

K = 1/2 (r + k) , Q = l/2(r - k) , R = l/2( gi + q 2 ) , 

C = R-^K-m B , B = Q .^ K , (5.7) 

and it satisfies the requirements of parity conservation and gauge invariance: 

M^k" = r»M^ = . (5.8) 

To calculate the matrix elements (5.2) and (5.5) in the DSB, we introduce two OVBs a^ and 
bjt (A = 0, 1, 2, 3) using the 4-momenta p\, p 2 , k and q 1 , q 2 , k: 

a o = P+/\Jp+ , a 3 = p-l\j-p 2 - , a 2 = [a ■ a 3 } x k/p , oi = [o • a 3 } x a 2 , (5.9) 

P± = P2 ± Pi , o±<5 = ai ± iSa 2 , 5 = ±1 , a 2 k = , a\ = a\ = a\ = — a 2 , = — 1 . 
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b = q + /Jqi , 63 = q-f^e , h = [b • h]*k/p , h = [b ■ b 3 ] x b 2 , (5.10) 

q± — q 2 ± qi , b± 5 > =b 1 ±i5b 2 , 6 — ±1 , b 2 k — , b\ — b\ — b\ — -b\ = -1 . 
where p and p are determined from the normalization conditions. Then the electron and proton 



j j 



operators P 21 ' and Q 21 ' [see (2.11) and (2.12)] will have the form 

P 2 { = 1/4 (m+pi) as a a* s , P 2 1 ' = <V 2 (m + pi) a s a 3 , (5.11) 

Q%f = 1/4 (M + gi ) by b % , Q/' 5 ' = S'/2 (M + gi ) b s > h , (5.12) 

while the matrix elements (5.2) and (5.5) in the case of various combinations of electron and 
proton spin states reduce to a product of traces: 

Mx = 1 Tr( P± S > 5 Q» ) Tr( Qff > S ' T ,(q 2 ) ) , (5.13) 



Mo 



\ Tr( P 2 f V ) Tr( Q^' ' 5 'm^ ) . (5.14) 






In the unpolarized case it is most efficient to use the calculation of the matrix elements in the 
DSB in conjunction with the standard approach [10]. The calculations performed by the first 
[i.e., using (5.13) and (5.14)] and second methods give identical result. Nevertheless, the second 
method, which will also be discussed below, is preferable, because it gives results considerably 
more quickly. To find the probability for the process (5.1) it is sufficient to calculate only the 
matrix elements of the electron and proton currents: 

(J^ S ), = u ±5 (P2hXPi) 5 = Tr( P 2 ^\ ) , (5.15) 

(J p ±5 '< 5 '), = u ±S \q 2 )T,{qV\qi) = Tr( Q 2 f '^(g 2 ) ) , (5.16) 

and also the quantity 

x ±6>,s> = u ±& \q 2 )M^u 5 \ qi ) = Tr( Qf? ' S ' M^e" ) . (5.17) 

The calculations give [11,36,47] 

(J 5 /), = 2m(a ) M , (J~ s '% = -25y4a s ), , (5.18) 

(J d p '> s '), = 2g e M(b ), , (j/' y ) M = -28'y'_g m {b 5 ,)» , (5.19) 



where y_ = J —p 2 _/2 , y_ = J—q 2 _/2, and g e and g m are just the proton electric and magnetic 
form factors [10]: 

q 2 

9e = fi + ^v^p h, 9m = /1 + V-p h ■ (5-20) 

Therefore, in the DSB the matrix elements of the proton current corresponding to transitions 
without spin flip are expressed in terms of the electric form factor g e , and the interaction with 
spin flip is expressed in terms of the magnetic form factor g m . 
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After the matrix elements of the proton current (5.16) are determined, the calculation of the 
contribution of the two Bethe-Heitler graphs reduces to the calculation of VCS on the electron 

[47,65]: 

\M, | =jMri{ J P _ 2pik z + e ^ J p )u( Pl )\ . (5.21) 

Denoting the result of averaging and summing the expression \M 1 ' | 2 over the polarizations of 
the initial and final particles by Y ee , we obtain [47,65]: 

F ee = 1/4 J2 Tr i (h + m) Qf># (px + m) Q? & }/Q 4 , (5-22) 

<5'e 



where Q^ 5 ' S — (Qe) (Jp S ' S )n * s the operator in parentheses between the electron bispinors u{p 



2 



_±<5',<5' 



and u(pi) in Eq. (5.21), and Q e = 70 (Qf 5 ' s ) + 7o- Owing to the factorization of the electric 
and magnetic form factors g e and g m in (5.19), the Bethe-Heitler term in the cross section for the 
reaction ep — > ep~y Y ee (5.22) will contain only the squares of the Sachs form factors (see Refs. 
36,47,65,67, and 68). 

Similarly, the calculation of the contribution of the graph in Fig. lc reduces to the calculation 
of quasireal Compton scattering on the proton. Using the expressions for the electron current 
(5.18), we have 

I M^ 5 | 2 = 1 I u(q 2 ) Qt 5 ' S u( qi ) | 2 , (5.23) 

where Q^ 6 ' 6 = (J^ s,s )^ M^ v e u . Denoting the result of averaging and summing Eq. (5.23) over the 
polarizations of the initial and final particles by Y pp , we obtain [65] 



Y pp = 1/4 ]T Tr{ (q 2 + M) Q^ 5 (q, + M) Q p }/r 4 , (5.24) 



— ±5 ' s -,„ 

where Q p = 7 (Q p ) + 'y°- Finally, to calculate the interference term in the case of unpolarized 

particles 

Y ep = 1/4 J2 2Re M i M 2 ( 5 - 25 ) 

<5,<5',e 

we shall use the matrix elements of the proton current (5.19) and also the 4- vectors X^ s ,s (5.17), 
which have the form [65] 

X^ ' = -25 y_b x k I ~^2~ T 2 + D2 T 4 + id y + y_— — — 2 — T 6 I e , 



2 



K* = 2 U (^ (r, + f^T 2 j + ^ (T, + ^T 4 ) ) + (5.26) 



-H <TpL v " — , " " T s \eT, 

where y+ = J q\/2 = My/1 — t , r = q 2 /AM 2 and z/i = kq + /2M 2 . As a result, for the matrix 
element M 2 (5.5) we find 

M 2 = w(p 2 )X ±,5 '' <5 'w(p 1 )/r 2 , (5.27) 
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and Eq. (5.25) reduces to calculation of the trace [65]: 



.±S'S' 



Y ep = 1/4 J^ 2Re {Tr ((p 2 + m)Qf' 5 '(p 1 + m)X ' )}/q 2 /r 2 , 



(5.21 



8,8' 



where X ±5 '> 5 ' = Y X± 5 '> 5 ' and X ' = (X± 5 '' 5 ')*^. The interference term Y ep (5.28) is a linear 
combination of the proton electric and magnetic form factors, because the operators Qf s ' s are 
expressed linearly in terms of the matrix elements of the proton current: Qf s ,s = (Q e )^(J p 5 ,s ) IX , 
[see Eqs. (5.3) and (5.19)]. 

Therefore, the problem of finding the probability for the reaction ep — > epj in this approach 
has been reduced to calculation of the traces (5.22), (5.24), and (5.28), which was done by means 
of the program REDUCE. For the differential cross section we then obtained [65]: 



da 



a 3 | rp 1 2 ^{j) X _|_ q 1 — p 2 — q 2 — k) d 3 p 2 d 3 q 2 d 3 k 

2p 2 o 2g 2 o 2w 



2n 2 Mp 1 q 1 ) 2 - m 2 M 2 



T \ 2 = 1/4 ]T | M fi \ 2 = F ee + Y ep + Y pp , 

pol 

8M 2 



m 



Yi 

2 



Y ee = ^ ( 9 2 e Y + T ^ Y U 

Ai A2 
A2 Ai 



m 2 g 2 / 1 



r 2 q 2 



2M 2 (l-r) V A 2 

Ai ^ A2 

A2 Ai 



Ai A2/ 2A1A2 

fPiQ+ P2q+\ 2 _ r {{pi q+f + (p 2 g+) 2 ) 
Ai ) (1-r) 



2 



m 2 g 2 / 1 



1 



A1A2 

r 2 q 2 



+ 



m 



piq^ 



2M 2 (l-r) V A 2 



m m 



P2q+ 
Ai 

2\2 



Ai A2/ 2A1A2 

r {{Piq+f + (P2q+) 2 ) 



2 1- — ) + ! m z ( — - t- 



1-T 



1 



A1A2 



K 



32M 3 



ep 



+r# 



r 2 q 2 (Av 2 — z/ 2 
1/1M 



Ax A 



fi-e^e 



Ai A- 



j/i 7i + T 2 



1/2 M: 



v x M 



1 -r 



1-r / V 1-r 

ite(j/iT 2 + 1/2T4) + 4MRe{ Zl T 2 + z 2 T 4 + ^ 3 T 6 ) 



Ta 



(5.29) 
(5.30) 
(5.31) 



(5.32) 



(5.33) 



(5.34) 



Y 



pp 



(aia 3 + 1/3) [(1 - r)|T\| 2 + 2z/ 1 Mi?e(T 1 T 2 *) + M> 2 - z/ 2 )|T, 



+ (a 2 + 1/3) [(1 - r)]T 3 | 2 + 2v 1 MRe(T 3 T:) + M 2 (v{ - u 2 )\T 4 



+ (a 2 o;3 + a 2 + 2z/ 3 )r 



|T 5 | 2 M 2 lml9 M 16M 4 



M 4 v% 



+ 



|T fi | 2 



«3 



(5.35) 
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For the invariant variables in Eqs. (5.30)-(5.35) used in determining the Bethe-Heitler term 
(Y ee ), the interference term (Y ep ), and the term corresponding to proton bremsstrahlung (Y pp ), we 
used the notation adopted in Ref. 64: 

yi = 2a 1 [a 1 a 3 (u 2 u 5 - v x v^) + 2v\ + u 2 u 3 ], v x = kq + /2M 2 , v 2 = -kq_/2M 2 , 

y 2 = 2a 2 (v 2 v 5 - 1/11/4) - «iZ/ 2 2 , 1/3 = r 2 /4M 2 , 1/4 = kq+/AM 2 , 1/5 = p+g+/4M 2 , 

y 3 = -(4z/ 3 /z/ 2 ) [aia 3 (^i^(^2 + ^3) - 2i/ 4 (i/ii/ 4 - 1^5)) + v A (Av 2 A - i/f)] , 

«i = ^5 + ^i^(2z/ 3 + z/ 2 )/z/f , a 3 = z/f/(z/f + (z/ 2 + vz){vl - v\)) , 

a 2 = m 2 /M 2 - 1/3 + M 6 /^ 2 [-(^i^4 + z^ 5 ) 2 + 4z/ 3 (z/ 2 - 1/11/41/5) - 4z/ 3 z/ 4 2 (z/ 2 + ^3)] , 

D 2 = M 6 [v\ + (z/ 2 + z/ 3 )(z/ 2 - y\)) = M 6 u 2 2 /a 3 , Ax = Pl k , A 2 = p 2 k , 

z-i = i>iv A a\a 3 , z 2 = v 2 v±a 2 , z 3 = l/Aai (2v 2 (2a 2 + v 2 + 1^3) + 4z/ 2 - z/ 2 ) . 

We note that the expression obtained for the differential cross section (5.29) coincides, apart 
from the definition of the initial quantities (the tensor M^), with the result obtained in Ref. 64, 
if in the latter f\ and f 2 are expressed in terms of g e and g m . Nevertheless, the Bethe-Heitler 
term Y ee and the interference term Y ep have a more compact form, owing to the factorization of 
the electric and magnetic form factors. 

Let us consider the effects due to contribution of all three graphs to the cross section for the 
reaction (5.1) in the selected kinematics when the initial proton is at test [qi = (M, 0)] and the 
electron beam energy is E e = 200 MeV. Performing the required integration over the phase space 
in the rest frame of the initial proton, we obtain [65] : 

o?ijj 2 \ q \T \ 2 
da = t — — dE pk dVt q2 dtt^ , (5.36) 

lbii z M I pi I {p 2 kj 

where <if2 7 and dQ q2 are the elements of the photon and proton solid angles, and E p k is the kinetic 
energy of the recoil proton. 

Let calculate the differential cross section (5.36) numerically in the region 5 < E p k < 35 MeV 
with the sum and the difference of the electric (a p ) and magnetic (/3 P ) polarizabilities equal to 
®p + P P = 14 and a p — f3 p = 10 (in units of 10~ 4 /m 3 ) [58-60]. We assume that the reaction 
kinematics is planar, and that the photon emission and proton scattering angles are $ 7 = 135° 
and d p = —20.5°. (All angles are measured from the direction of motion of the primary electron 
beam). The calculation [65] show that in the entire range of proton kinetic energy considered, 
5 < E p f, < 35 MeV, for the selected angles i? 7 = 135° and -d p = —20.5° the electron scattering 
angle $ e and the 4-momentum transfer \r\ = J— (p 2 — p\) 2 are bounded by the values |i? e | < 6.4° 



and \r\ < 7.3 MeV, with the minimum value of |r| corresponding to forward electron scattering. 

The results of numerical calculations of the differential cross section (5.36), da / dE p k/ dVL q2 / dVL 1 
in the kinematics described above are shown graphically in Fig. 2. We see that in the angular 
range studied the cross section for the reaction ep — > ep~f has a sharp peak consisting of two 
maxima. This peak originates from the factor 1/r 4 in Eq. (5.35) for Y pp . The two maxima 
have a kinematical origin and arise from the interference of two pole graphs corresponding to 
quasireal Compton scattering. The cross section (5.36) has a strong angular dependence, which, 
in particular, causes the two maxima to disappear when the proton (or photon) emission angle is 
changed by only one a degree (i.e., for -& p = —19.5°), so that we have an ordinary peak at E p k = 25 
MeV. 
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Figure 2: Differential cross section (5.36) for reaction ep — > ep~f in the kinematics where proton 
bremsstrahlung dominates, see comments in the text. Proton scattering and photon emission 
angles are i? p = —20.5° (solid line), i? p = —20.0° (dashed line), -& p = —19.5° (dot-dashed line), and 



ih 



135°. 



ep I / ^pp 

—20°) give results 
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The differential cross section (5.36), shown by the graphs in Fig. 2, is the sum of the Bethe- 
Heitler (a ee ), the interference (a ep ), and the proton (<%>) terms [see (5.30)], where the symbol (a) 
denotes a cross section of the form (5.36) with | T | 2 replaced by Y ee , Y ep , and Y pp , respectively. 
Numerical calculations shown that in the entire range of proton kinetic energy studied, 5 < 
E p k < 35 MeV, the ratios of the Bethe-Heitler term a ee and the interference term a ep to the term 
corresponding to proton emission o vv are bounded by the values (T ee /a pp < 0.02 and |ov„|/<t„„ < 
0.05. The calculations carried out for another set of angles (i? 7 = 135° and d p = 
which are only insignificantly different: <7 ee /<jp P < 0.05 and |<7 ep |/cr pp < 0.075. 
are much smaller than unity, the main requirement (see Ref. 60) for separation of the background, 
which is mainly electron bremsstrahlung, is satisfied. 

To explain the sensitivity of the reaction ep — > ep~y to the proton polarizability we performed 
numerical calculations of the cross section (5.36) for the same set of angles (i? 7 = 135° and 
dp = —20°) for fixed sum of the electric and magnetic polarizabilities a p + (3 P = 14 but different 
values of the difference: (a) a p — (3 P = 10 and (b) a p — f3 p = 6. It turned out that the cross section 
(5.36) is about 8% larger for the smaller difference of polarizabilities. Therefore, in this kinematics 
the cross section for the reaction ep — > ep'j is quite sensitive to the proton polarizability [65]. 



6 Emission of a linearly polarized photon by an electron 
in the reaction ep — » ep^ 

Let us consider the emission of a linearly polarized photon by an electron in the reaction ep — > 
ep'-f, taking into account the proton recoil and form factors. Our study will be limited to the 
contribution of the two Bethe-Heitler graphs (a) and (b) in Fig. 1, which corresponds to the 
matrix element (5.2). The contribution of the graph with VCS on a proton can be neglegted when 
the initial electrons have ultrarelativistic energies, and the photon and final electron are scattered 



27 



at small forward angles (i? 7 ~ m/E e , $ e ~ m/E e , m/E e <1). 

We are interested in these effects for the following reasons. First, even though the Bethe- 
Heitler process has been studied earlier in the case of the emission of linearly polarized photons 
[69,70] and is widely used to obtain them at accelerators [71], up to now the proton recoil and 
form factors have not been accurately taken into account (in contrast to the unpolarized case). 
Second, as was shown in Ref. 72, the inclusion of these factors in the case of unpolarized photons 
leads to a strong change of the differential cross section for the Bethe-Heitler process. Since the 
polarization characteristic of the scattered radiation are expressed in terms of the differential cross 
section for the emission of an unpolarized photon (see below), it is clear that inclusion of the recoil 
and form factors is essential. 

The covariant expression for the differential cross section for the Bethe-Heitler process (in the 
Born approximation) taking into account the proton recoil and form factors in the case of emission 
of a linearly polarized photon has been obtained by us in Ref. 73. It has the form 

a 3 | T e | 2 5 4 (p 1 + qi~p 2 -q 2 -k) d 3 p 2 d 3 q 2 d 3 k 
d&BH = , , (o.IJ 

2tt 2 yWi) 2 - m2M2 2p20 2q20 2uJ 

AM 2 
\T e \ 2 =^ r {g 2 e Y I z + TglY i n, (6.2) 

A 2 Ai i — t A\A 2 

Yl ? = -2-±-± + -!- ^±£ + (q 2 + 4m 2 ) (ea) 2 - 4 (eA) 2 , (6.4) 

\ 2 Ai i — t \i\ 2 

Pi P2 A , %Pl)P\ (b Pl)P2 , a _v 

a = — - — , A = b + . (6.5) 

A\ A 2 A\ A 2 

All the quantities entering into (6.1)-(6.5) are defined in the preceding section. Thus, the differ- 
ential cross section for the Bethe-Heitler process in the case of emission of a linearly polarized 
photon doBH (6-1) naturally splits into the sum of two terms containing only the squares of the 
Sachs form factors and corresponding to the contribution of transition without (~ g 2 Y} e ) and 
with (~ r g^Yjf) proton spin flip. 

Let us discuss the properties of the 4- vector a, which is well known from the theory of emission 
of long- wavelength photons [10], and the 4- vector A. They both satisfy a condition which follows 
naturally from the requirement of gauge invariance: 

a k = A k = , 

and, in addition, they are spacelike vectors: a 2 < and A 2 < 0. This is easily verifed by using 
the 4- momentum conservation law and the explicit form of a 2 and A 2 : 

a2=m2 {{-i) 2+ jk' (6 - 6) 

,2 _ i , m2 (<i+Pi i+pA , T q+Pi-q+Pi / fi7 x 

A _1 + 4M 2 (l-r) \T\ 2 XT) + T^F X 1 \ 2 ■ {bJ) 
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We note that the 4- vector A was first introduced in Ref. 73. 

Using the electron 4-momenta p\ and p 2 and the photon 4-momenta k, we construct the 4- 
vectors of the photon linear polarization en and e± (enk = e±k = e\\e± = 0): 

(p 2 k)p 1 - (pik)p 2 [pi ■p 2 } x k 

ej| = ; , e± = ; , (b.8) 

P P 

where p is determined from the normalization conditions: e| = e\ — — 1. Then the degree of 
photon linear polarization will be given by the following expressions [73]: 

_ | T ± | 2 - | T|| | 2 _ A x 
7 ~ I T ± I 2 + I 71, I 2 ~ A ' (6 ' 9) 



where 



16 M 2 

A = — (gj A n + r g 2 m A 12 ) , (6.10) 

Q ?M"2 

A 2 = ^ r (g 2 e Y 1 + rg 2 m Y 2 ), (6.11) 

4 = A 2 + rM 2 a 2 +2(eJ ) 2 , (6.12) 

Ai 2 = -A 2 + r M 2 a 2 - 2(e±b ) 2 + m 2 a 2 , (6.13) 

±{SDf 
M 2 (l-r)a 2 A?A| 

W<x(Pi) M (P2)^(gi) p (g ; 
r, = 2 - 2i _ 2£ _ _L_ i^±l! _ 2 rM 2 a 2 - 2 A 2 , (6.16) 



Mo) 2 = - ^; V 7,^2 , (6-14) 



M 2 (l -r)a 2 A?A1 

5D = 1/2 e, I/pa (p 1 r(p 2 nq 1 y(q 2 r , (6.15) 

Ai A 2 t (kq^ y -' 

X 2 Ai 1 — t A1A2 

F 2 = -2 - ^1 - ^ + -^- &±^ -2rM 2 a 2 + 2A 2 -2m 2 a 2 . (6.17) 

A 2 Ai 1 — T A1A2 

It is easy to check that A 2 (6.11) coincides with the expression for Y ee (5.31) determining the 
Bethe-Heitler cross section in the case of unpolarized particles: A 2 = Y ee , and also that Y\ = Yj 
and Y 2 = Y n [see (5.32) and (5.33)]. 

Therefore, owing to the factorization of the squares of the form factors g e and g m and also the 
use of the 4- vectors a and A (6.5), the differential cross section for the Bethe-Heitler process in 
both the cases of linearly polarized photon (6.2) and unpolarized photon (6.11), (5.31), can be 
written in a rather compact form. 

Let us integrate Eq. (6.1) over d 3 q 2 and dp 2 Q in the rest frame of the initial proton, qi = (M, 0). 
As a result, we find: 

dcr B H _ a 3 uj \p 2 \ \T\ 2 
dojdQjdQe ~ (2vr) 2 \p x \ q 4 ' l ' j 

\T\ 2 = g 2 e Y I e + rg 2 m Y I f. (6.19) 

Let us consider the limit of the cross section (6.18) when the proton is a point (structureless) 
particle with infinite mass, i.e., we assume that g e = g m = 1 and q 2 = (M,q) ~ (M, 0), where 
q = p\ — p 2 — k is the momentum transferred to the proton. In this limit [M — ► 00), E^ v = 
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q 2 /2M -> 0, q/2M -> 0, and 6 = (1, g/2M) ~ (1,0). We choose the Coulomb gauge for the 
photon polarization vectors: e = (0, e), as a result of which we find 

Pie p 2 e pie p 2 e , , 2 2 2 

e6 = 0, ea = — , eA = p 20 p w — , r{q + k) = u q . 

Ai A 2 Ai A 2 

Using these expressions to take an limit in (6.19), we obtain: 

Ai A 2 uo 2 q 
A 2 Ai AiA 2 



T\ 2 = 2-^-^-^ + q 2 (ea) 2 + 4 (eA) 2 , (6.20) 



or, in expanded form, 



b. _ bi - u2q2 ( 4 2 2 ) C Ei^ x 2 

A 2 Ai AiA 2 V Ai 



+ ( 4p?o + g 2 ) (^) 2 - 2 ( 4 PloP20 + g 2 ) ^-^ . (6.21) 

The expressions (6.18) and (6.21) for the differential cross section for the Bethe-Heitler process 
do bh I 'duj j 'dVl 1 / 'dVt e in the limit where the proton is an infinitely heavy, structureless particle 
coincide with the analogous expressions of Ref . [69] . 

7 Virtual-photon polarization in the reaction 

ep — » epj (ep — » eX) 

The reaction ep — > ep~y and VCS on a proton have recently become interesting not only at low and 
intermediate energies [60], but also at high electron energies and 4- momenta transferred to the 
proton [63,74-77]. The VCS process offers greater possibilities for studying hadronic structure than 
the RCS process, because in it the energy and three-momentum transferred to the target can be 
varied independently. These attractive properties of VCS have led to the suggestion that it be used 
for experimental study of the nucleon structure [74,75] and have made it necessary to perform a 
thorough theoretical study of the reaction ep — > epr) (including the use of the noncovariant method 
of calculating helicity amplitudes; (see Refs. 63, 76 and 77 and references therein)). To calculate 
VCS on a proton, it is necessary to know the hadron (W^ u ) and lepton (L^) tensors [63,78]: 

L,u, = JpJZ, J M = u(p 2 )i tl u(p l ) , (7.1) 

where u(pi) are electron bispinors, u(pi)u(pi) = 2m, and m is the electron mass (i = 1,2). The 
interpretation of the results is considerably simplified if the tensor L^ v is expressed in terms of 
the longitudinal and transverse polarization vectors of the virtual photon. The corresponding 
expressions can be found in Refs. 63 and 78. However, they have two defects: (1) the electron 
mass is neglegted, which is of course justified at ultrarelativistic electron energies and large squared 
4-momentum of the virtual photon; (2) they have a noncovariant form. A lepton tensor free of 
these defects was constructed in Ref. 79. 

Let us consider the question of the polarization state of a virtual 7 with 4-momentum r — p\ — p 2 
which is exchanged between the electron and proton in the reaction ep — > epj (see Fig. lc). Using 
the vectors of the orthonormal basis a^ (5.9) (A = (0, 1, 2, 3)): 

a = P+/\Ip 2 + , a 3 = p-/\j-p 2 - , a 2 = [a • a 3 ] x qi/p , a x = [a • a 3 } x a 2 , (7.2) 
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p±= p 2 ±pi , a 2 qi = , a\ = a 2 , = a 2 = -a 2 , = -1 , (7.3) 

which satisfy the completeness relation 

Oo ' °o — a i ' a i ~~ a 2 ' 0,2 — a 3 • 03 = g , (7-4) 

we construct the 4-vectors of the longitudinal (63) and transverse (ei, e 2 ) polarization of a virtual 
photon with 4-momentum r (Ref. 79): 

[o -Oi]?i [a -a 3 ] x q 1 (1 + a 3 ■ a 3 )q 1 

ei = . = , e 2 = a 2 = , e 3 = = , (7.5) 



>3?i) 2 + g? p VM1) 2 + gf 

where 

2 / n2 2j9ip 2 -j9 1 g 1 -p 2 gi -M 2 ((pip 2 ) 2 -m 4 ) -m 2 ((pi<?i) 2 + (p 2 <?i) 2 ) 

P = (oigi) = j; r^ i . 

(P1P2) ~ m 

It is easily verifed that the 4-vectors e^ (i = 1, 2, 3) are orthogonal to each other (e^e, = 0, i ^ j), 
and also that e^r = e^ = and e 2 = e 2 , = — e 2 = — 1. The 4-vectors e* (7.5) are not changed 
when the auxiliary 4- vector g x is replaced by 91 + pi — p 2 = g 2 + k [because p\ — p 2 = r = —2ya 3 , 
where y = \ / —r 2 /2, and because the vectors a a (7.2) are orthogonal]. For this reason, study of 
the virtual-photon polarization vectors e* (7.5) in the rest frame of the incident proton or in the 
cm. frame of the final proton and photon is equivalent and leads to the usual expressions. Here 
we shall restrict ourselves to the rest frame of the incident proton \q\ = (M, 0,0,0)], where the 
4-vectors e^ have the form: 

e 1 = (0,1,0,0), e 2 = (0,0,1,0), e 3 = -^=(\ f |,r„n 3 ) . (7.6) 

Here n 3 is a unit vector directed along r {n 2 = 1), and tq is the time component of the 4- vector 
r = (r ,r). 

The four mutually orthogonal vectors ei,e 2 ,e 3 , and 03 also satisfy the completeness relation: 

e 3 -e 3 -e 1 -e 1 -e 2 -e 2 -a 3 -a 3 = g , (7.7) 

which allows a^ and a,\ to be expressed in terms of e\ and e 3 : 

a\ = ae 3 — fit\ , a = f3e 3 — ae\ , j3 2 = 1 + a 2 , (7.8) 

aiQi n a oQi / 7n x 

a = e 3 a± = a^e\ = . , p = e±ai = e 3 ao = . . (7.9) 



(a 3 gi) 2 + q\ v (a 3 gi) 2 + g 2 

In the DSB (4) the matrix elements of the electron current have the form of (5.18): 

(J 5 /), = 2m(a )„ (J- S > s )„ = -25y (a 5 )„ (7.10) 

where a±s — 01 ± i5a 2 , S — ±1. Let us write them in terms of the 4-vectors e; (7.5) (Ref. 79): 

(J^V = 2m (/3e 3 - aej, , (J^' 5 )^ = -2^ (ae 3 - /3ei + %8e 2 )^ . (7.11) 
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Therefore, for transition without electron spin flip (Jf' 5 ) the virtual-photon polarization vector 
is a superposition of the longitudinal (fie$) and transverse linear (— ae\) polarizations, while for 
transition with spin flip («/ ( T <5 ' <5 ) it is a superposition of the longitudinal (aes) and transverse 
elliptical [eg = (0, e^) = —f3e\ + i5e 2 ] polarizations. Here the state of a photon with elliptical 
polarization vector e& = (0, eg) will have degree of linear polarization (equal to the ratio of the 
difference and sum of the squared semiaxes [57]) [79]: 

„-£=!--*- (712) 

Inverting this relation, we obtain: 

p^l + K^ ^ Q 2 = J^7_ (713) 

Now we find the squared moduli of the vectors e$ and a$: 

| e s \ 2 = 1 + {3 2 = -?— , | a 5 \ 2 = (l + (3 2 ) (1 + k l ) , (7.14) 

1 K"y 



'0 -»2 _ '0 

- r 2) ' e 3 ~~ f- r 2^ 



«£ — ft 7 e 3 — k 7 / 2 \ ' e 3 — ^ „2"\ ' ('-15) 

We introduce the normalized vectors eg ' and as ': 



e , =7rT^ = V^ e5 ' |e ^ =l - (7 ' l6) 



a& = /1 , g2 = V 2 ^'^l = 1 + K 7 e 3 = ! + K i , ( 7 - 17 ) 

Therefore, the elliptical-polarization vector e^ of a virtual photon can be normalized to unity 
(\e$ '\ 2 = 1), but the presence of a longitudinal polarization makes this normalization impossible 
for the total vector a$ ' simultaneously. The quantity k l (7.15) corresponding to the inequality 
\a§ '\ 2 = 1 + k,l 7^ 1 has the meaning of the degree of longitudinal polarization of a virtual photon 
emitted in a transition with electron spin flip. In the ultrarelativistic limit, when the electron 
mass can be neglected, the quantities k 7 and Kl will be interpreted as the total degrees of linear 
and longitudinal polarization of the virtual photon. In this (massless) case we have: 

Mi) 2 + q\ = -M 2 ^ , (a iqi ) 2 = M 2 ctgH/2 , (7.18) 

K- l = l-2 r —tg 2 §/2, (7.19) 

where i? is the angle between the vectors f>\ and j9 2 . Equation (7.19) for k 7 coincides with the 
result of Ref. 78. 

The vector as ' (7.17) can also be written as 



as =V^i n 3-\/ — 7T^ e i+™\ — tt~^ e 2 , (7.20) 
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which makes it easy to construct the polarization density matrix for a virtual photon in the 
massless limit (both in the polarized case, which for massless particles is helical polarization, and 
in the unpolarized case; see Ref. 78). 

To obtain the complete expression for k 1 and Kl arising from the contributions of the matrix 
elements both without and with spin flip, we construct the lepton tensor averaged over electron 
spin states. Using the matrix elements (7.10) and (7.11), this can be done fairly simply [79]: 

L^ v = Am 2 {a Q ) IJ ,(a ) v + % 2 ((ai)^(oi) v + (a 2 )^(o 2 ) v ) • (7.21) 

Using the completeness condition (7.3) and gauge invariance, the tensor L^ v can be written as 

L^ = Ax 2 (a ) M (ao)^ - % 2 SV , (7.22) 

where x 2 = m 2 +y 2 . The tensor L^ v (7.22) can be used to reduce the calculation of the contribution 
of graphs with VCS on a proton to the cross section for the reaction ep — > epj to calculation of 
the trace of a product of tensors: 

Y pp = V W^ , W^ = VX , V„ = u(q 2 ) M^e" u( qi ) 1 . (7.23) 

Let us express the tensor L^ u (7.21) in the terms of the virtual-photon polarization vectors e, (7.5). 
As a result, it naturally breaks up into the sum of three terms corresponding to the contributions 
of transverse (Lt) and longitudinal (Ll) states and their interference (Llt) [79]: 

L = Ay 2 (L T + L L + L LT ) , (7.24) 

L T = ei • ei ((3 2 + a 2 m 2 /y 2 ) + e 2 • e 2 , (7.25) 

L L = e 3 • e 3 (a 2 + (3 2 m 2 /y 2 ) , (7.26) 

Llt = ~ (ei • e 3 + e 3 • d) a/3 (1 + m 2 /y 2 ) . (7.27) 
Then the total degree of linear polarization of the virtual photon will be given by 

= (3 2 + a 2 m 2 /y 2 - 1 = a 2 

7 (3 2 + a 2 m 2 /y 2 + 1 (3 2 + 1 - 2m 2 /a; 2 ' l ' J 

Since a and /3 are the same in Eqs. (7.12) and (7.28) [see (7.9)], the inclusion of the electron mass 
in the ultrarelativistic limit will lead only to a slight increase of k 7 [79]: 



Inverting the relation in (7.28), we find 



-L t K~, o ^o o / o ^fv^ 



m 



2 



^ + a 2 m 2 /y 2 = — -i , « 2 + (3 2 m 2 /y 2 = —^ + ^ . (7.30) 

l — k 1 l — k^ y 

We can separate the completely polarized and unpolarized parts in the transverse tensor Lt (7.25): 

2 
Lt = ei-ei (,3 2 +a 2 m 2 /t/ 2 -l)+e r ei+e 2 -e 2 = ( k 7 e r ei+(l-K 7 ) (e r ei+e 2 -e 2 )/2 ) . (7.31) 

1 K~f 
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Therefore, the virtual-photon polarization density matrix pij is obtained from the tensor L„ (7.24) 
just as in the massless case (see Ref. 78): 

Pij = (1 - «y L l3 /8y 2 . (7.32) 

For the degree of longitudinal polarization of the virtual photon we then obtain: 

1- 




(7.33) 

The expressions (7.28) and (7.33) for k' and k' l with m = obviously become k 1 and /t^ of (7.12) 
and (7.15). 

We conclude by noting that the region of applicability of the tensor L^ v (7.24) is not limited to 
only VCS on a proton. Since in fixed-target experimets the charged- lept on scattering at available 
energies is mainly determined by virtual photon exchange, the tensor L^ (7.24) can also be used 
to study deep-inelastic electron scattering (e ± p — »■ e ± X), and muon scattering (fi^p — > fJ^X), 
where inclusion of the mass is more important. 

8 Compton back-scattering of the photons of a circularly 
polarized laser wave on a beam of ultrarelativistic, 
longitudinally polarized electrons 

It was shown in Refs. 80 and 81 that, using existing (SLC) and planned (VLEPP) accelerators 
with colliding e + e~ beams, it is possible to obtain colliding 7e and 77 beams of roughly the same 
energy and luminosity as the original e + e~ beams. It have been suggested that the intense beams 
of hard 7 rays needed for this be obtained from the Compton back-scattering (CBS) of a powerful 
laser flash focused on the electron beam [82]. For a sufficiently powerful flash in the conversion 
region [81], processes with simultaneous absorption of several laser photons from the wave become 
important: 

e~ + n 70 — ► e~ + 7 , n > 1 , (8.1) 

7 + s 7o ->■ e + + e" , s > 1 . (8.2) 

The first of these nonlinear processes leads to broadering of the spectrum of high-energy photons 
[83], and the second effectively lowers the e + e~-pair production threshold [84]. 

The process (8.1) and (8.2) were studied systematically in Ref. 85. In Ref. 16 they were studied 
from the view-point of providing sources of polarized 7 and e + e~ beams. The phenomena arising 
in collisions of polarized electrons with the photons of a circularly polarized electromagnetic wave 
were analyzed in Ref. 86. Nonlinear effects were studied not only for £ 2 < 1, but also for £ 2 > 1. 
Here x 2 is the wave intensity parameter: 

e=nJ^), (8.3) 

where n 7 is the photon density in the wave and uo is the photon energy, a is the fine structure 
constant, m is the electron mass. The emission spectra at high intensities (£ 2 > 1) were first 
calculated numerically in Ref. 83, but the particle polarization was not taken into account. 
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Recently at the SLAC accelerator a series of experiments [87] are being performed for £ ~ 1 to 
verify nonlinear QED. This has become possible owing to the use of supershort, strongly focused 
laser pulses. The region of nonlinear effects for £ 2 > 1 is very important here, and it is of great 
interest because emission processes due to simultaneous absorption of a large number of photons 
from the wave become important, and the probabilities for these processes are essentially nonlinear 
functions of the field strength. 

As a rule, in the literature the laser wave is described as the field of a planar electromagnetic 
wave [85,86]. The applicability of this model in strong fields has been studied in Ref. 88. 

According to Ref. 10, the S'-matrix element for the transition of an electron from the state 
i>p — V ;5 (P) S ) to the state ip p i = i[) ±s (p',s') ,(5 = ±1), with the emission of a photon of 4- 
momentum k! = (u/, k') and circular-polarization vector ey is given by 

S fi = -ie f^ p , e* x , ip p exp(ik'x) (2cu')~ 1/2 d A x , (8.4) 

where ip p and ip p , are the exact wave functions of electrons in the field of a circularly polarized 
electromagnetic wave, corresponding to the vector potential 

A — ai cos(A;:r) + \a,2 sin(kx), A = ±1. (8-5) 

Here k is the wave vector, k 2 = 0, a\k = 0,2k = a\d2 — 0, a\ — a 2 , — a 2 , and A, A are helicities 
of a laser and emission photons. The explicit form of the matrix elements (8.4) in the DSB was 
obtained in Refs. 38 and 86: 



ie (4 



7T 



,1/2 



% - ^ rojr'Lw* £ M ±8 ( 2 -) 4 5 > k + Q-J-V), (8-6) 

[Zuj zq zq ) 1 n=1 



(n) _ 1 w ( _ x)n c{_ 2(1 ~ U/Un ) 

2 l ' M y/w' - 1 



M_S = - A' (-A)» £ - ^^ ( Jn _ x + Jn+l] 



+ ~r^ ( ^H - 8A> / ) J n+ J , (8.7) 

2(U+1) VVW'+1 y/w'-l) n+AA J' V > 



(n) 1 , / / , ,„ >. / U /_ U \ ( U + 2 \VV' - 1 



^■' = - 5 V(-A).{W - 1- - — VwTT-^' 



where 



x (v TTF^- 1 + Jn+l) " 2(J+\)Vw- 1 Jn+AA 'j ' (8 ' 8) 

^ ?TZ , / t 771 2 /2 2 2 /1 1 /-2\ 

q=p+^—k,q=p+ 1 ^- T k ) q =q = m^ = m (1 + £ ) , 
(1 + AA') T (1-AA') T , ,, , 

«/n+AA' = ?; Jn+1 H ~ -Ai-l , uk + q = k + q , 



2n£ / w / w 



Vl+£2 V M « ^ M r 
35 



(n) 



(n) 



Here M & s and M_ s s are the emission amplitudes of the n-th harmonic corresponding to tran- 
sitions without and with electron spin flip, q and q' are the electron quasimomentum 4-vectors, 
1 — (Qo,Q), q' = (q'o, q'), and J n are the n-th-order Bessel function of argument z n . It is eas- 
ily verifed that the amplitudes M±$ s have the following kinematical features. For u = u n and 
n > 1 they vanish [M±g S (u = u n ) = 0]. The reason for this behavior of the amplitudes will 



be explained below. Knowledge of the diagonal amplitudes (8.7) and (8.8) allows transformation 
to the helicity amplitudes (see Ref. 86). As a result, we obtain the following expressions for 
the differential cross section of the hard photon emission by an electron in the field of circularly 
polarized electromagnetic wave [86]: 



ira 



da c __ 

du xm 2 ^ 2 (u 



7 , ( Fi n + AA e F 2n + AA F 3n + A e A F, 



n=l 



An 



J.10) 




fS.ll^ 



It, 



'in 



M + l 



;-4 .11 + e (2 + m) {Jl_x + J 2 n+1 - 2 J; 



here x = 2kp/m 2 , A e is the helicity of the electron, A e = ±1. The expression inside the summation 
in (8.10) determines the emission probability of the n-th harmonic when the polarization states 
of the laser and the emitted photons and also the initial state of the electron are helicity states. 
For £ 2 = Eq. (8.10) coincides with the result of Ref. 89. 

Using (8.10), the degree of circular polarization of a photon in the final state A/ is defined as 



A, 



£ (A F 3n + A e F 4n ) / £ ( F ln + AA e F, 



2n 



(8.12) 



n=l 



n=l 



For £ 2 < 1 only the first few harmonics dominate in the cross section (8.10) for the process (8.1). 
We expand the expressions (8.11) in the parameter A = £ 2 /(l + £ 2 ), expanding only the Bessel 
functions and using the exact expressions for the u n . As a result, for the first three harmonics we 
have [86]: 
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for the third harmonic. 

The inclusion of the third harmonic, whose probability is proportional to A 2 , leads to the 
appearance of terms containing A 2 in Eqs. (8.13) and (8.14). This is the main difference between 
the result obtained in Ref. 86 for the emission probability of the first two harmonics and the 
analogous expressions from Refs. 16 and 85. 

Let us consider the case of a head-on collision of ultrarelativistic electrons with the photons of 
a laser wave. To obtain the energy distribution of the produced photons da c /dy, where y = uj'/E, 
and E is the electron energy, in (8.10) we must make the replacement u — ► m/(1 — y) [85]. Here 
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variation of the variable u in the range < u < u n correspond to variation of y in the range 
y ■ < y < y n , where 

nx u n nx 2kp AuiE 

1 + £ 1 + m„ nx + 1 + £ z m z m z 

Comparing the maximum possible energy of photons produced in ordinary Compton scattering 
(n = 1, £ 2 = 0) with the energy calculated with inclusion of nonlinear effects (£ 2 7^ 0), we see that 
photons of the first harmonic (n = 1) have lower maximum possible energy. However, the energy 
of 7 quanta emitted in the absorption of several photons (n > 1 + £ 2 ) is greater than that available 
in ordinary Compton scattering. Making the replacement: u —* y/{l — y) in (8.10) and (8.11), we 
obtain the distribution in the energy of the hard 7 quanta y = uj'/E [86]: 

da no 2 °° 

-T- = —^H E ( F i« + AA - F 2n + AA' F 3n + A e A' F 4n ) , (8.16) 

ay xm 2 £, 2 ^ 

F ln = -AJ 2 n + e (l- y + 1 ±-\ ( Jl_ x + J 2 +1 - 2J 2 n ) , 




'2n = S "ityt 1 - ^" 7^ -7- Wn-1 - ^n+lJ > 



^3n=^ U-yt- \ l ~ Z H V y J n-l~ J n+l) 1 



'In 



(8.17) 



V / T+F V "' n Vn\ Vn) (l-l/) 2 ' 

Let us now turn to the more detailed analysis of the influence of nonlinear effects on this 
process. We shall start from the following initial. We take a head-on collision to be one in which 
the electrons have energy E = 50 and 300 GeV, and uo = 1.17 eV (a neodymium laser). We shall 
use the expansions (8.13)-(8.15) for numerical calculations of the energy spectra (l/W) dW/dy 
(where W = J2n=T W n is the total emission probability) and the degree of circular polarization 
Xf of an emitted photon for £ 2 < 1. For ( 2 > 1 we shall use the exact expressions (8.16) and 
(8.17). In this case n max is determined from the conditions for the series (8.16) to converge. 

The results of numerical calculations of the energy spectra for various polarizations of the 
initial electrons (A e ) and laser photon (A) are shown by the graphs in Figs. 3a, 3b, and 3c 
for £ 2 = 0.3, 1, and 3, respectively. We see from these figures that the inclusion of nonlinear 
effects leads to a significant difference between the calculated spectra and the spectra of ordinary 
Compton scattering. First, the simultaneous absorption of several photons from the wave leads to 
broadering of the hard-7 spectrum and the appearance of additional peaks corresponding to the 
emission of higher-order harmonics. For a given electron energy this broadering is larger, the larger 
the wave intensity. For example, for E = 50 GeV and £ 2 = 0.3 the spectrum is bounded above by 
the value y ~ 0.67, while for £ 2 = 1 it practically vanishes at y ~ 0.8, even though an insignificant 
fraction of the photons can carry off up to 97% of the electron energy. Second, the effective increase 

38 



of the electron mass [85] m 2 



m; 



m 



£ 2 ) leads to compression of the spectra at smaller 



values of y, because for each n the spectrum is bounded above by the value y n = nx/(l + nx + £ 2 ) 
and not by nx/(l + nx). The increase of the electron energy decreases the relative compression 
of the first harmonic (see Fig. 3a). For relatively low intensity of the laser wave (£ 2 = 0.3) the 
main contribution to the emission comes from photons of the first harmonic, and the yield of 
photons from higher harmonics is insignificant. At intermediate intensity (£ 2 = 1) the broadering 
of the spectrum due to nonlinear effects is accompanied by an increase of the probability, and the 
yield of harder photons becomes important. Finally, at high intensities (£ 2 = 3), as seen from 
Fig. 3c, emission owing to nonlinear multiphoton absorption processes becomes comparable to 
one-photon emission and even begins to dominate (at E = 50 GeV). Therefore, emission of the 
first harmonic dominates in the CBS spectra in the field of a circularly polarized electromagnetic 
wave at £ 2 = 0.3, while at £ 2 = 3 the emission in mainly due to higher harmonics, i.e., the emission 
of a hard photon by an electron essentially becomes nonlinear [86] . 
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Figure 3: CBS spectra correspond for the following values of the intensity parameter £ 2 : (a)— > 
0.3; (6) — > 1; (c) — > 3. The dashed lines correspond to ordinary Compton scattering (£ 2 = 0). The 
lines 1, 2, and 3 correspond to the following choice of helicities of the electron and laser photon: 
1 _+ A e = 0, A = 1; 2 -> A e = 1, A = -1; 3 -+ A e = 1, A = 1. 

To study the polarization effects at each value of the energy E, we calculated the energy spectra 
for the following polarization states of the electron and laser photon: 



1 



A, = 0, A = 1: 2 -» A, = 1, A = -1: 3 -+ X P = 1 , A = 1. 



These correspond to lines 1, 2 and 3, respectively, in Fig. 3. Everything said above about the 
behavior of the energy spectra pertained to these three lines. Regarding their relative location, 
from Fig. 3 we see that the most intense spectra correspond to the case where the electron and 
laser photon spins are parallel (AA e = — 1), while the least intense ones correspond to antiparallel 
spins (AA e = 1), as in the case of ordinary Compton back-scattering (see Ref. 89). 

We also note that the difference between the spectra calculated for the three polarization cases 
considered is very large at small values of the intensity parameter (£ 2 = 0.3), but insignificant at 
£ 2 = 3 (E = 50 GeV). It again arises only in connection with increasing electron energy (see Fig. 
3c for E = 300 GeV). 
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Figure 4: Energy dependence of the degree of circular polarization of high-energy photon, calcu- 
lated at £ 2 = 0.3 for the following polarization states of the colliding particles: (a) A e = 0, A = 
1; (b) X e — 1, A = — 1; (c) A e = 1, A = 1. The dashed lines correspond to ordinary Compton 
scattering. 
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Figure 5: Energy dependence of the degree of circular polarization of high-energy photon, calcu- 
lated at £ 2 = 1 for the following polarization states of the colliding particles: (a) A e = 0, A = 
1; (b) X e — 1, A = — 1; (c) A e = 1, A = 1. The solid lines correspond to electron energy E = 50 
GeV, and the dashed lines to E = 300 Gev. 



Let us consider the energy dependence of the degree of circular polarization of a hard 7 ray, 
shown by the graphs in Figs. 4 and 5. For this we first note that the above-mentioned kinematical 
features of the behavior of the amplitudes Mj^j in (8.7) and (8.8) has a spin origin [86]. In 
fact, the equation u = u n correspond to photon emission in the direction of motion of the initial 
electron beam. In the case of absorption of n photons (n > 1) from the wave and exact backward 
scattering of the hard photon, the total helicity of the e + rr/o and e + 7 system before and after 



the interaction is not conserved. It is this which causes all the amplitudes M± s 



(«), 



U = U r 



for n > 1 



and also M, 



W, 



8,5 



U 



Ui)to vanish. The requirement of helicity conservation also leads to Aj 



-A 



for ordinary Compton scattering at the edge of the spectrum [86]. 
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As can be seen from Figs. 4 and 5, the inclusion of nonlinear effects (£ 2 7^ 0) decreases the 
degree of circular polarization at the first peak. The contribution of higher harmonics leads to the 
appearance of additional peaks, and at the edge of the spectrum (for n = n max ) we have A/ = —A, 
as in the case of ordinary scattering. However, it should be noted that the yield of these photons 
is insignificant, since the spectra are practically broken off at y <C y nmax - The situation regarding 
AA e = — 1 is the most favorable in this respect, as there is a large range of hard 7 energies in 
which the degree of circular polarization | A/ | is very close to unity. 

9 e + e~-pair production by a hard photon in a 
collision with photons of a laser wave 

In Ref. 84 it was shown that a hard photon obtained in the reaction (8.1) can create e + e~ pairs in 
a collision with photons of the same laser beam. The threshold for this reaction (8.2) at s = 1 is 
very high. The lowest energy of the Compton photon (s = 1) in the process (8.2) for a neodymium 
laser with uo = 1.17 eV is u = m 2 /ujQ = 223 GeV. In fact, e + e~ pairs will be created in large 
numbers and at significantly lower energies owing to collisions of the hard photon 7 with several 
laser photons 70 simultaneously [84]. Observation of the process (8.2) is particularly interesting 
for verifying QED in a new parameter region. At the same time, it is an important source of 
background for 7e and 77 collisions, and a possible method of dealing with it is described in [84]. 

Like (8.1), the reaction (8.2) is an interaction of electrons and photons with the field of an 
electromagnetic wave which is nonlinear in the field strength. It is easily checked that the inclusion 
of the influence of the nonlinear effects in (8.1) on the process (8.2) also leads to a significant 
lowering of the e + e~-pair production threshold and to an increase in the number of pairs [90]. 

The maximum energy of a Compton photon 7 resulting from the absorption from the wave of 
n laser photons of energy c^o by an electron of energy E is 

nx AujqE 

u n = — E , x = — . 9.1 

1 + nx m z 

The threshold value of the 7 energy for the process (8.2) is given by 

(k + sk ) 2 = Am 2 , (9.2) 

where k and ko are the 4-momenta of the photons 7 and 70. The corresponding threshold values 
of the energy of the electrons in the accelerator beam E ns for e + e - -pair production owing to 
absorption of n photons from the wave and collisions with s laser photons are determined from 
(9.1) and (9.2): 

E ns = ^—(l + {l + s/n) l l 2 ). (9.3) 

For n = 1 we obtain Eq. (7) of Ref. 84. Using (9.3), we can calculate the values of Ei s and E 2s 
for c^o = 1-17 eV and 1 < s < 6. The results (in GeV) are given in Table 1: 

These results clearly show that the broadering of the hard-7 spectrum due to nonlinear effects 
also leads to lowering of the e + e~-pair production threshold. 
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Table 1: Threshold values of the electron energy in the accelerators beam E ns (in GeV) for e + e 
pair production at various n and s in the case of neodymium laser 



s 


1 


2 


3 


4 


5 


6 


E\s 


269 


153 


112 


90 


77 


68 


E 2 s 


248 


135 


96 


76 


64 


56 



The matrix elements M±^„ = MJ.^ and the differential probability for the process (8.2) in 
the field of a circularly polarized electromagnetic wave are given by [90]: 



dW (s) 



e 2 m 2 

Atiuj 



M^' \ 2 6\sk + k - q - q>) 
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, d 3 qd 3 q' 
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(9.4) 
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Here ho, A and k, A' are the 4-momenta and helicities of the laser and hard photons, // is the pro- 
jection of the positron spin on the axis (1.9), q and q are the positron and electron quasimomenta, 
so is the threshold value for the number of absorbed photons, J s = J s (z) is the Bessel function of 
argument z, and £ 2 is the wave intensity parameter (8.3). 

The total probability for pair production by a photon in the process (8.2) per unit volume and 
unit time is given by [90]: 



W 



am 

4w 



2 oo 



]T / (F 0s + XX'F 2s + fiXG 0s + fiX'G 2s ) 

s>s i 



du 



uJu(u 



(9.5) 
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*o. = Jl + ? (2u - 1) (- J 2 S + ( JU + J s 2 +1 )/2 
F 2s = e (2u - l)(2u/u. - WU - J 2 +1 )/2 , 

G Qs = ij+^j- £ 2 u/u s (u s - 1) (J s 2 _ x 



Js+l ) ) 



G 2s = ^-{u s J 2 s + e[usJ 2 s + u{u - 1){JU + J 2 s+l ) ]} , 
V± = l/v/(w'±l)/2. 

The total number of e + e~-pairs N e + e - created by a hard photon is obtained by summing over 
the energy of the Compton photons [84] : 



N P 



^if/^-Wa)^^ 



where iV 7 is the total number of hard photons, <J C (E) and da c / 'duo are the total and differential 
cross sections for CBS, W^(u,uq,^) is the probability for pair production by a hard photon per 
unit time in the process (8.2), and r is the duration of the laser flash. The results of numerical 
calculations of log(N e + e -/N e ) as a function of the beam electron energy E for various energies of 
the laser flash A, wave polarizations A, helicity of the initial electron beam A e , and spin projection 
/i on the C3 axis (1.9) for positrons are shown by the graphs in Fig. 6 [90]. 
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1 50 E, GeV 




Figure 6: Dependence of the number of e + e~-pairs created by a hard Compton photon on the 
electron beam energy. The lines 1, 2, and 3 correspond to the following choice of helicities A and 
A e of the initial particles in the reaction (8.1): (1) AA e = 0, (2) AA e = — 1, (3) AA e = 1. The 
solid lines correspond to n — 2, and the dashed lines to n — 1 in (8.1) Figure (a) corresponds to 
positron spin projection /j, = — 1, and (b) to n — +1. 

The lines 1, 2, and 3 correspond to the following choice of helicities: (1) AA e = 0; (2) AA e = — 1; 
and (3) AA e = 1. The solid lines correspond to n — 2, and the dashed lines to n — 1 in the process 
(8.1). It follows from Fig. 6 that nonlinear effects in CBS lead to a significant increase in the 
number of e + e~ pairs created by a hard photon at current accelerators energies. 
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Conclusion 

The goal of the present review was to explain what the DSB is, what new contribution it makes 
to the description of particle spin properties and also to the development of the covariant method 
to calculate matrix elements in the Bogush-Fedorov approach, and how this method is related to 
others like that of the CALCUL group. 

We have seen that the DSB plays a key role among all the other methods in that in it the 
Loretz little group common to particles with 4-momenta p\ (before the interaction) and p^ (after 
the interaction) is realized. The DSB allows the description of the spin states of systems consisting 
of two particles (even when they have different masses) by means of the spin projections on a single 
common direction. The coincidence of the Loretz little groups causes the particles before and after 
the interaction to have a common set of spin operators which commute with each other, and this 
allows the covariant separation of the interactions with and without change of the spin states of 
the particles involved in the reaction, so that the dynamics of the spin interaction can be traced. 
Thanks to the coincidence of the spin operators and also the fact that Wigner rotations are singled 
out, the mathematical structure of the diagonal amplitudes is maximally simplified. 

To calculate matrix elements in the covariant Bogush-Fedorov approach, it is necessary to know 
the projection operators of the particle states, the operator for the transition from the initial to 
the final state (and its inverse), and also the raising and lowering spin operators in the case of 
spin flip transitions. In the review we have developed this covariant approach by using the DSB. 
We have constructed the operators u s (pi) u ±5 (p3) used to calculate the diagonal amplitudes in 
the case of transitions without and with spin flip. They are valid in both the massive and the 
massless cases. We have obtained three equivalent representations for them which have a compact 
form. We have also studied the transition to the massless case, in which the DSB coincides up to 
a sign with the helicity basis. 

In the CALCUL method the fermion must be massless. The key feature of that method is the 
very convenient choice of photon polarization vectors, in which the momenta of the fermions from 
which the photons are emitted are used. This ensures gauge invariance and simplifies the structure 
of the amplitudes, so that they can ultimately be calculated. In the CALCUL method the mass 
can be taken into account only in the ultrarelativistic case and only in the form of awkward mass 
corrections. As a rule, generalizations of this method to the massive case require the introduction 
of auxiliary vectors unrelated to the kinematics of the problem, which are therefore inconvenient 
to work with. Nevertheless, this method contains the attractive idea of constructing the photon 
polarization vectors in terms of the 4-momenta of the particles participating in the reaction. 
This allows a decrease in the number of various scalar products in the final expressions for the 
amplitudes and thereby simplifies the calculations. Therefore, giving up on the generality of 
the treatment makes the solution of the problem more efficient. This is even more true with 
regard to the method developed for calculating diagonal amplitudes, because the construction of 
the mathematical formalism for them involves only the 4-momenta of the particles participating 
in the reaction. In the DSB this is sufficient, thanks to the use of the ideas of the covariant 
Bogush-Fedorov approach. 

Let us briefly list the main results of our calculations of several specific QED processes using 
the method developed for calculating matrix elements in the DSB. 

We have shown that in the ultrarelativistic (massless) limit, the differential cross sections for 
Moller and Bhabha bremsstrahlung (e ± e~ — > e ± e~ r y) in the case where not only the initial e^ and 
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e _ , but also the photon are helically polarized can be presented as the product of two factors, 
one of which is universal and coincides with that obtained earlier by the CALCUL group when 
polarization is absent. 

The helicity amplitudes of the three-photon annihilation of a free pair e + e~ — »■ 37 have been 
calculated along with the orthopositronium annihilation amplitudes corresponding to total spin 
projection 0, ±1. The differential cross sections taking into account the polarizations of the various 
particles were obtained. The annihilation probability was calculated in the case where one 7 is 
linearly polarized and the other two are unpolarized. The expression obtained for the degree of 
photon linear polarization coincides with the results of other authors. 

A compact expression was obtained for the differential cross section of the Bethe-Heitler emis- 
sion of a linearly polarized photon by an electron, taking into account the proton recoil and form 
factors, thanks to the factorization of the squared electric and magnetic form factors of the proton. 
In the limit where the proton is a point particle of infinite mass, this expression becomes the usual 
one. 

We have studied the reaction ep — > epy, taking into account the proton polarizability in 
the kinematics corresponding to electron scattering at small angles and photon scattering at 
fairly large angles, where proton bremsstrahlung dominates. The results of numerical calculations 
performed in the rest frame of the initial proton at electron beam energy E e = 200 MeV in 
the chosen kinematics show that the conditions needed to isolate the subprocess '-/p — * 7p from 
the reaction ep — > epy are satisfied, because the relative contribution of the Bethe-Heitler and 
interference terms to the reaction cross section is less than 10 %, and the cross section for the 
reaction ep — > ep~y is quite sensitive to the proton polarizability. 

A covariant expression has been obtained for the lepton tensor in which the contribution of 
states with transverse and longitudinal polarization of the virtual photon is isolated. It has been 
shown that inclusion of the lepton mass tends to increase the degree of linear polarization of the 
virtual photon. 

We have studied nonlinear effects in Compton back-scattering of photons by an intense circu- 
larly polarized laser wave focused on a beam of longitudinally polarized ultrarelativistic electrons 
(e + n'-fo — * e + 7). We have found that at high intensities the emission of a hard photon is 
essentially nonlinear, and the effect of the polarizations is markedly diminished. 

We have shown that the broadering of the spectrum in nonlinear Compton back-scattering 
tends to lower the e + e~-pair production threshold and increase the number of pairs in collisions 
of a hard Compton photon with several laser photons simultaneously (7 + rr/o — > e + + e~). 

Thus, the absence of difficulties associated with inclusion of mass and calculation of spin flip 
amplitudes, and also the elegance of the results obtained, demonstrate the clear superiority of 
developed method to calculate matrix elements. 
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